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1. Introduction 

The story of the Bohr collective Hamiltonian seems to go back to 1879, long before 
the discovery of atomic nuclei, when Lord Rayleigh [1] showed that the time dependent 
coefficients, a\^(t), in equation 



of the surface of an incompressible liquid drop in the spherical coordinates R, 6, play 
the role of normal modes of small oscillations of the surface around a spherical shape. 
The next milestone was the liquid drop model of an atomic nucleus [2] — an object 
known already at that time. It allowed Fliigge [3] to apply the Rayleigh normal modes 
to a classical description of low-energy excitations of spherical nuclei. The model of 
Rayleigh and Fliigge was quantized by Aage Bohr [4] who formulated in this way the 
quantum model of surface vibrations of spherical nuclei. It was Aage Bohr also, who 
introduced the concept of the intrinsic frame of reference for the quadrupole nuclear 
surface and replaced the variables a 2(1 by the Euler angles and parameters /3 and 7 
(nowadays often called the Bohr deformation parameters) which describe the surface 
in the intrinsic system. Afterwards, Bohr and Mottelson [5] generalized the model to 
vibrations and rotations of deformed nuclei. A generalization of the Bohr Hamiltonian 
to describe large-amplitude collective quadrupole excitations of any even-even nuclei 
was proposed by Belyaev [6] and Kumar and Baranger [7]. This general form of the 
Bohr collective Hamiltonian is sketched in Volume 2 of the Bohr-Mottelson book j8j. In 
the meantime several specific forms of the collective Hamiltonian intended to describe 
the collective excitations in nuclei of different types were proposed [9l ITO l [TTl [T2l LT~3l [14] . 
All that belongs to the past. During the last decades of the twentieth century and 
in the beginning of the twenty first century, a great progress in the development of 
microscopic many-body theories of nuclear systems took place. However, it is still not 
possible for the self-consistent Hartree-Fock-Bogolyubov approach, employing single- 
particle degrees of freedom, to directly describe the collective phenomena in nuclei. One 
of the methods which can be used in order obtain such a description is the Adiabatic 
Time Dependent Hartree-Fock-Bogolyubov theory (ATDHFB) [HI [15], [16] , which in the 
case of the quadrupole coordinates leads to the Bohr Hamiltonian. On the other hand, 
the Random Phase Approximation (RPA) (cf [T7|, [18] for its application to nuclear 
physics) is not able to explain a large amplitude collective motion because it assumes 
harmonicity of the excitations. The other method which can extend a microscopic 
approach to the collective phenomena is the Generator Coordinate Method (GCM) 
[j~9] . In its Gaussian Overlap Approximation (GOA) [201 EE] the GCM also yields the 
Bohr Hamiltonian [22, 23j. The integral Hill- Wheeler equation of the GCM without 
this approximation for all five quadrupole generator coordinates is considerably more 
difficult to solve [24]. 




(1.1) 
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It is not an intention of the present article to review all applications of the Bohr 
Hamiltonian used to describe the nuclear quadrupole collective excitations which have 
been done hitherto. Neither is the full list of publications on this subject compiled here. 
Rather the review focuses on the Bohr Hamiltonian itself. Section [2] contains a detailed 
discussion of the kinematics of the quadrupole degrees of freedom. Tensor properties 
of the coordinates and momenta themselves are a subject, which goes far beyond the 
Bohr Hamiltonian model. The discussion is based on the concept of the isotropic tensor 
fields presented in |Appendix A[ The tensor algebra and analysis allow us to investigate 
in section [3] the possible most general form of the Hamiltonian and other observables. 
Moreover, the tensor structure of the collective wave function with a definite spin is 
presented. In section [3] we also recall some analytically solvable examples of the Bohr 
Hamiltonian and we briefly review different methods of solving the eigenvalue equation 
of the general Hamiltonian. A specific version of the collective model is determined by 
several scalar functions which define the Hamiltonian and all other observables. Within 
the phenomenological approach to the collective model the parameters of all these scalar 
functions are fixed using experimental data on properties of the collective excitations. 
Obviously, there is another possibility, that all the necessary tensor fields, such as the 
collective potential, the inertial functions, the moments of charge distribution and the 
gyromagnetic tensors, are derived from an underlying more fundamental theory. Among 
others, the classical or semi-classical theories or models, such as different versions of 
the liquid drop model (e.g. [HI [25J [26]) or the Thomas-Fermi model |27l 28] can be 
used and because of such a possible classical background the Bohr Hamiltonian is 
usually qualified as the geometrical model. However, we focus here on the quantum 
methods of derivation of the Bohr Hamiltonian from a microscopic many-body theory. 
In section [4] the ATDHFB theory, which leads to the classical (at the intermediate stage) 
collective Hamiltonian, and the Generator Coordinate Method, which gives the quantum 
Hamiltonian directly, are recalled and their application to the collective quadrupole 
Bohr Hamiltonian is presented. The section ends with an example of calculations which 
illustrates the entire reviewed approach that leads from a microscopic theory through the 
Bohr Hamiltonian to the description of the collective quadrupole excitations in several 
even-even nuclei ( 178 Hf — 200 Hg). Some important properties of the ATDHFB inertial 



functions are discussed in Appendix B 



2. Quadrupole collective coordinates 

The fundamental assumption of the Bohr collective model applied to the description 
of the quadrupole collective states is that the dynamical variables (coordinates) form 

(E) 



a 



a real quadrupole electric, i.e. of the positive parity (cf |Appendix A.lj ), tensor 
(the superscript (E) will be omitted below). It is not necessary to give the tensor «2 
a geometrical interpretation. In particular, the interpretation of «2 coming from fll.ll) . 
even if often used, is not needed. Thus, referring to the model as to the geometrical 
model is, in general, unjustified. Indeed, at the present time the Bohr model with 
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some special forms of the collective Hamiltonian is formulated and treated simply as an 
algebraic collective model (cf e.g. [29 l l30l IBT] ). 



2.1. The laboratory coordinates 

Let the covariant components a 2fl (/•* = — 2, . . . , 2) of a tensor a 2 in the laboratory 
system Ui a b form a set of five independent dynamical variables. Although in general, 
it is convenient to use complex coordinates, there are, in fact, five real coordinates 
02(U (l 1 — 0, 1, 2) and b 2fJi (fi = 1, 2) which are respectively the real and imaginary parts 
of a 2fJi (see (I A. 71 ) for definition). The volume element in the space of these coordinates 
is (cf flSH) 

2 2 

dfi= n d ^n d& 

The covariant momentum tensor and the contravariant momentum Hermitian adjoint to 
it (cf ( 1A.22I) and (1A.23fl ) are n 2fJi = —ihd/da 2 and 7rL = —ihd/da 2 ^, respectively. The 
angular momentum vector connected with rotations in the physical three-dimensional 
space, expressed by the five components of the quadrupole tensor, reads (cf (IA.26I) ) 

Liv = L { $ =h/T0[a 2 x 7t 2 ] lv . (2.2) 



2.2. The intrinsic coordinates 

Apart from, perhaps, the case of the five-dimensional harmonic oscillator Hamiltonian 
[4], the laboratory coordinates are not used in practice. Usually, one introduces the 
system of principal axes of the tensor a 2 , Ui n , called the intrinsic system. The orientation 
of U; n with respect to Ui a b is given by the three Euler angles, (p, ip. For given 
laboratory components a 2fl the Euler angles are defined through the three implicit 
relations: 

1 £ t D% :1 (<p,0,1>)a 2lt = O, (2.3a) 



^D%{ V ,^)a 2il = 2^D» 2 (<p,#,il>)a 2li (2.36) 

which imply vanishing of the following real and imaginary parts of the intrinsic 
components a^fl of a 2 : 

off = = fo| n) = 0. (2.4) 
The two remaining real parts of the intrinsic components are 



(in) _ 

ho = a o 



£ tf^W, (2.5a) 



4 n) = a 2 = -L {D%(^ 0, il>) + D 2 ;_ 2 (if, 0, ^j)) a 2 ,. (2.56) 
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These coordinates are usually parametrized by two parameters, (3 and 7 (0 < (3 < 00, 
~n < 7 < 7r) defined by the relations 

a = /^cos7 a 2 = /?sin7. (2.6) 

The parameters /3 and 7, commonly known as the Bohr deformation parameters, play 
the role of polar coordinates in plane (a a 2 ) (see figured]). In equations (12.3a!) . (12.361) . 
(12.5 al) and (12.561) the Wigner functions with v 7^ appear in the following combinations 
denoted as 

tf&faW) = (DlM$^) ± (-l)^(^^,^)) • (2.7) 
Relations (I2.3aj) . (12.361) . (12.5 al) and (12.561) define the transformation 

a2 M (a , a 2 , y?, tf, -0) = D^fo, tf, ^)a + D 2 ^\ip, tf, V0 a 2 (2.8) 

from the intrinsic coordinates ao, a 2 , y?, V" to the laboratory ones — a 2 ^ (/i = 
—2, . . . , 2). The Jacobian of this transformation is, up to a constant coefficient, equal 
to (cf [32]) 

D(a 2 - 2 , . . .,0*22) 



D(a , a 2 , (p,fi,ip) 



\X(a ,a 2 )\ sintf (2.9a) 



D(a ,a 2 ,v2,^,^) 

X(a ,a 2 ) = a 2 (3aQ — a 2 )a 2 (v^3a + a 2 )(v / 3a-o — a 2) = /3 3 sin 37. (2.96) 

The transformation (12.61) to the polar coordinates (3 and 7 adds the factor (3 to the 
Jacobian (12.961 ) since da da 2 = (3d(3&). For values of the intrinsic coordinates for which 
the Jacobian is not equal to zero the transformation is invertible, meaning that the 
intrinsic frame can be defined. The inverse transformation is given by implicit functions 
of the Euler angles. The Jacobian vanishes on the three straight lines in the (a a 2 ) plane 
specified by the conditions 

a x = or 7 = 0, ±7r a y = or 7 = ir/3, — 27r/3 a z = or 7 = 27r/3, — n/3 (2.10) 
where 

a x = -- (V3a + a 2 ) = (3 shit,, a y = - (VSa - a 2 ) = /3 sin7 y 

a z = «2 = /3sin7 2 and 7^ = 7-2^/3 7^ = 7 + 2^/3 7 2 = 7. (2.11) 

The lines intersect at the point ao = a 2 = {(3 = 0) and divide the plane into six sectors 
as is shown in figure [TJ At the boundaries of these sectors the intrinsic frame cannot be 
uniquely defined. The reason is that the tensor a 2 is then axially symmetric (cf ( 1A.2I) ). 

As will become clear, the appearance of the three lines, instead of only one, on 
which the Jacobian vanishes, reflects the fact that equations (12.3al) and ( 12.361) define, in 
general, three mutually perpendicular axes in the three dimensional space but not their 
ordering (say x, y, z) nor their arrows. This means that the axes are defined up to a 
discrete group of 48 transformations O*,, k — 0, . . . , 47, changing the names and arrows 
of the axes. The group is the symmetry group of a cube and is denoted by Oh (cf [33]). 

Traditionally, one generates all the 48 transformations Ok G Oh by three so called 
Bohr's rotations, Ri, R 2 , R3, and the inversion P (other sets of generators can also 
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Figure 1. Three straight lines: a x — 0, a y = 0, a z — (or six rays: 7 = 0, 7 = ±7r/3, 
7 = ±27r/3, 7 = ±7r) going through the origin of the coordinate system and dividing 
the (0002) plane into the six sectors. Each point of the polar coordinates (/3o,7o) 
located in the sector < 7 < it/3 has five twin points in the remaining five sectors. 
All these six points of polar coordinates (/3o, ±7o), (A), ±27r/3 ± 70), marked with full 
circles, are equivalent to each other. This means that any intrinsic component of a 
tensor field at one of these points is a linear combination with constant coefficients of 
the intrinsic components of the same field at any other of the remaining five points. It 
is therefore sufficient to consider tensor fields in one sector only. Usually, the shadowed 
sector, < 7 < it/3, is considered. 

Table 1. Generators of the Oh group of transformations of the coordinate system. 
The index k numbers the transformations. The angles w^, Qk, Cfc are the Euler angles 
defining the orientation of the coordinate system after the transformation Ofe. The 
coordinates Xk, Uk, Zk define the position with respect to the transformed system. 



Transformation The Euler angles Coordinates 



fc 


Description 


o fc 




Qk 


ft 




Vk 







Identity 


E 











X 


y 


z 


1 


Bohr's 


Ri 


7T 


7T 





X 


-y 


—z 


2 


rotat- 


R 2 








tt/2 


I) 


—X 


z 


3 


ions 


R.3 





tt/2 


tt/2 


I) 


z 


X 


24 


Inversion 


P 











—x 


-y 


—z 



be chosen). Definitions of the Bohr's rotations are given in table [TJ The twenty four 
unimodular transformations (rotations) Ofc for k = 0, ... , 23, which are the products 
of the Bohr's rotations are listed in [35] (but the numbering in [35] is different from 
that used in this paper). The remaining twenty four transformations can be formed by 
composing the unimodular ones with the inversion. 
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Table 2. Transformation rules of a 2fl and b 2fi under Ri, R2 and R3. The real and 
imaginary parts of the a 2 components transformed by Rfc are denoted by a^} and b 2 k J . 
Decomposition into two irreducible representations of the group O, E and F 2 , is also 
shown. 

E F 2 






E 


a 2 o 


a22 


a 2 \ 


b 2 i 


022 


1 


Ri 


a-20 


0,22 


-021 


621 


-022 


2 


R 2 


a 2 o 


-a 22 


&21 


-&21 


-022 


3 


R 3 


— ^0,20 + ~^T a 22 


\/^3 1 
2~Q20 — 2° 22 


-021 


b 22 


-«21 


24 


P 


«20 


022 


a 2 i 


621 


b 22 



The transformation rules of the tensor a 2 under the Oh transformations are 
presented in table El Since the tensor a 2 has the positive parity, it is insensitive to 
the inversion of the coordinate system. For a 2 it is therefore sufficient to consider 
only the SO (3) group of rotations. The five-dimensional irreducible representation 
of SO (3) can be decomposed into two irreducible representations, E and F 2 , of the 
octahedral group O. These irreducible representations are identical with representations 
of Oh [33]. This is clear from table [21 where the transformation rules of the real and 
imaginary parts of a 2fl under the Bohr's rotations are shown. The bases of the irreducible 
representations E and F 2 are easily identified. Clearly, a 20 and a 22 form a basis of the 
two-dimensional irreducible representation E of O, whereas the components — a 2 i, — b 2 i, 
b 22 form a basis of the three-dimensional irreducible representation F 2 . This makes 
obvious that equalities (12.41 ) are invariant under the Oh transformations. 

When names and arrows of the intrinsic axes are changed, intrinsic coordinates cto, 
a 2 (or /3, 7), <f, ip of (|2.8p are transformed accordingly. Changes of the deformation 
parameters and of the orientation angles under the O transformations are given in table 
[U Values of the Euler angles (p 3 , $3, t/> 3 after the circular permutation R 3 of the intrinsic 
axes can be found from the following equations given in [34J: 

cos $3 = — sin $ cos ip sin $3 sin ip3 = cos $ cos (^3 — ip) = cosip sin ^3. (2.12) 

The Wigner functions of the Euler angles (pk,$k,?Pk giving the orientation of the 
rotated intrinsic system are related to those of the Euler angles <p, ip of the orientation 
of the initial intrinsic system by the following superposition formula |34j: 

Dl v {pk^k^k) = Y; D lAv^^) D l>^k,Q k ^k) (2.13) 

where Wk, Qk, * with k — 0, 1, 2, 3 are the Euler angles of the corresponding rotation. 

In order to present the contravariant and covariant momentum tensors as differential 
operators in the intrinsic variables, one should first calculate the derivatives of the 
intrinsic coordinates with respect to the components of the tensor a 2 using relations 
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Table 3. Transformations of the Bohr deformation parameters and of the Euler angles 
defining the orientation of the intrinsic axes with respect to the laboratory axes under 
the Bohr's rotations Rfc, k = 1,2,3, of the intrinsic coordinate system. The Euler 
angles ip3, $3, ip3 giving the orientation of the intrinsic system with axes permuted 
circularly can be found from equations l|2.12p . 

Bohr's rotation Deformations The Euler angles 

k Ok Wk Qk ft Pk lk p>k J?fc ipk 

0E000/3 7 (p <d *p 

lRl 7T 7T 7 If ±11 IT — d —1p 

2 R 2 tt/2 /3 -7 Lp <d V + tt/2 

3 R 3 tt/2 tt/2 7 -2tt/3 ip 3 d 3 ip 3 



(12.3d) . (12.351) and (I2.5al) . fl2.5 6j) . The calculation is a bit tedious. Relations (I2.3al) and 
( 12.351) . when differentiated with respect to c^, yield a system of three linear equations 
for the derivatives of the Euler angles from which the derivatives themselves can be 
calculated. Useful formulae for the derivatives of the Wigner functions with respect to 
the Euler angles can be found in [34]. By differentiating relations (12.5 al) and (12.561 ) with 
respect to a 2A1 one obtains the derivatives of a and a 2 . The complete calculation is 
performed in [32]. The final result for the momentum tensor reads: 

= [Dl (^ W)£- q + D% +) fa 0, ^) + " 2 £ ^ ! L z fa 0, V) 

* ^' - , — L v {<p,W) (2.14) 



where defined by ( 12.111) . and 

A / cos if) d d d 

L x (<P,ti,i/>) = - ihl — —7- + sin ^^ + cot0cosT/>— 
\ sini; oip ov oip 

f sin ib d d 3 \ 

L y (ip,$,ip) = -ihl- — -— + cos?/;— - cot$shi7/>— (2.15) 
\sinw dip ov dip J 

d 

L z {cp,$,iP) = ~ih— 

are the Cartesian intrinsic components of the angular momentum. Indeed, from (1A.26I) . 
(12. 8p and (12.141) one can derive the following relations 

Z$ = L M = =f-L (4 ± L,) (2.16) 

between the laboratory and the intrinsic spherical components of the angular momentum 
vector: 

K = T, D l^- (2-17) 

V 

It should be remembered that the spherical intrinsic components of angular momentum 
fulfill the commutation relations with the sign opposite to that of ( 1A.29I ). namely 

[LffM?] = +hV2(l^v\lK)Lt } (2.18) 
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and that these components commute with the laboratory ones 

[L^\L lv ] = 0. (2.19) 

Using the deformation parameters instead of the intrinsic components one should replace 
in (12.140 the derivatives with respect to ao and a 2 by the derivatives with respect to (3 
and 7: 

d d 1 d d d Id 

COS7— — sin 7— — — — — = sin 7— + COS7— — . (2.20) 



da d(3 (3d"f da 2 8(3 (3 8^ 

It should be noticed that both the laboratory and the intrinsic components of angular 
momentum do not commute with the Wigner functions and that the order of D^ v and 
L>x^ in (12.141) is essential (this is not so in (12.17ft — see [36]). However, it is sometimes 
convenient to exchange their order using the commutation relations. Doing so in the 
formulae for the Hermitian adjoint momentum one obtains: 

+ (£ + H£ - h - y ) ( D ^^r) - k iM) 

x y 

Using the deformation parameters (3 and 7 one should appropriately express the 
derivatives and take the dependence of a and a 2 on (3, 7 into account. The following 
relations are useful to this end: 

yl(±_l\ = 3sin2 ^ I (1 _ 1 _ 1^ = 3cos2 7 (2 22) 

2 \a y a x ) (3 sin 37 2 \ a z a x a y J (3 sin 37 
At the end of the discussion of the intrinsic coordinates, let them be functions of the 
parameter t (time in the case of classical motion). The laboratory coordinates become 
then the functions of t too. The derivatives of a 2 ^ with respect to t read: 

a 2 » = aoDl Q (<p, ip) + d 2 D^ 2 +) (p, <&, ip) 

-2\uo x a x D 2 X ] ((p, #, 1/)) - 2u; y a y Dl{ +) (<p, 0, V) + 2iu z a z D%- ) (<p, 0, V) (2.23) 
where 

uj x = •& sin ip — ip sin •& cos ip u y — i3 cos ip + (p sin 1? sin ip uj z = ip + p cos d (2.24) 

play the role of the intrinsic components of the angular velocity. The derivatives do and 
a 2 can be expressed by f3 and 7: 

d = (3 cos 7 — j3j sin 7 d 2 = /3sin7 + /37COS7. (2.25) 
£ 5. Functions of the coordinates 



The general rules presented in Appendix A for constructing isotropic tensor fields of a 



given tensor are applied here to the case of the quadrupole tensor a 2 . In this case one 
has the following five elementary tensors [37]: 
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(i) two quadrupole tensors 



■yJT/2[a 2 x a 2 \ 2ll oc e^(a 2 ) 



-2|t 1 
(2), 



«2 



(ii) two scalars 



-(2). 



[a 2 ■ a 2 ) oc e o{a 2 ) 

-(3), 



(7 = 

Xo = (a 2 ■ or 2 ) oc £$(a 2 ) 
(iii) one octupole tensor 

= V2[a 2 x a 2 ] 3fl oc ef^(a 2 ). 

The elementary tensors are related to each other by the single syzygy which 
[38]): 

IXs x X3\ 6M = ~y {2%o [a 2 X a 2 x a 2 ] 6M 

-3<r [« 2 x a 2 x a 2 ] 6M + [<r 2 x a 2 x t7 2 ] 6M } • 

Since the intrinsic components c4±i = ( see ( )2.3ah ). the only nonvanishing 
components of the elementary tensors are 



(2.26a) 
(2.266) 

(2.26c) 
(2.26 d) 

(2.26e) 
reads (cf 



a 



(in) 
20 



a = P cos 7 



(in) 
2±2 



^° 2 = ^ ,5sin7 



(T. 



(in) 
20 



2 2 «2 o 
a — a 2 = P COS ^7 



(in) 
2±2 



(in) 
X32 



(in) 
"X3-2 



T.V+o - v' 2 (•/()(•/ 2 

=a 2 (3aQ — a^) = —^/3 3 sin 37. 



--^=/? 2 sin 2 7 



^-v-u -2; ^ 
The scalars cr and x are the following functions of a and a 2 , or /3 and 7: 

cr = cto + a 2 = Xo — a o( a o — 3o3,) = (3 3 cos 37. 

Scalar functions depending on <7o and %o can be treated as functions of (3 and ( 
The fundamental tensors of a given even rank 2L (L 
elementary tensors a 2 and a 2 as follows: 



0, . . .) are built 



Tk2L = r 2 V 



a 2 x 

v 



■ ■ X a 2 X (T2 x 



X a 2 



2L-n 



n—L 



(2.27) 
intrinsic 

(2.28a) 
(2.286) 
(2.28c) 

(2.29) 

= cos 37. 
up from 

(2.30a) 



2L 



for k = 1, . . . , L + 1 (n = L, . . . , 2L), where A; stands for the tensor consecutive number 
and n — for its order in a 2 (n — L + k — 1). To construct the fundamental tensors of an 
odd rank 2L + 3 (the tensor with L = 1 does not exist) one aligns the even rank tensors 
(12.30aj) with the single elementary tensor \3 ( see (I2.27P ): 



_ (n+3) 
Tk2L+3 = T 2L+3 



a 2 x ■ - x 02 , x g 2 x - ■ x a 2- < xx3 

2L—n n—L 



(2.306) 



2L+3 



The ranges of / and n do not change. Obviously, for the fundamental tensors, only 
the intrinsic components with even magnetic numbers are nonvanishing. An analytic 
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formula for the intrinsic components of the fundamental tensors created by aligning 
several a 2 's alone is known [39], [40]. It reads: 

(in) 
21 2m 

where 

t2l2m{l) = 



Jm) 
' 1 21 2m 



a 2 x • • • x a 2 
v - ' 

l 



(2.31) 



(21 + 2m)\(2l - 2m)\ ( l \ ( r 
(47)! Km) V 6 



l—m 



\l—m I ~ 

cos 7J I —= sm 7 



V2 



x 2 Fi 



l — m l — m 

5 



-; m + 1; - tan 7 
3 



(2.32) 



2 ' 2 

with 2 F\ being the hypergeometric function [H]. Then, comparing (12.2861) with (12.28d) 
one finds that 



r, 



(in) 

1+121 2m 



02 X 



x cr 2 



(in) 



2i 2m 



^^(-27). 



(2.33) 



The fundamental tensor with a given even rank 2L and the number /c = n + 1 (n = 
0, . . . , L) is created by aligning the tensors ( 12.3ip and (12.331) with ranks 21 = 2L — 2n 
and 21 = 2n, respectively. To form the fundamental tensor of an odd rank, the third one, 
namely X3 °f (I2.28cl) should be aligned. For further use, it is convenient to introduce 
the 'dimensionless' intrinsic components of the fundamental tensors which depend only 
on 7: 



tn+12L 2K (7) = 
tn+1 2L+3 2K (7) 



(in) 



oL+n 'n+l2L2K 



[*2£-2n(7) X <2n(-27)] 



2L2K 



L+n+3 T n+12L+3 2K 



(in) 



Wl2i(7) X jpX3 



(2.34a) 
(2.346) 



2L+3 2K 

An arbitrary isotropic tensor field of the quadrupole tensor a 2 has the following 
structure: 

ki 

Ti M {ptz) = ^2 C) r fc/Af(«2, o- 2 , xs) (2-35) 



fc=l 



where 



1/2 + 1 
(J-3)/2 + l 



for even / 
for odd J 



(2.36) 



and fk are arbitrary scalar functions. There is no vector (tensor of rank 1) field of a 2 . 
In the special case of the quadrupole field (7 = 2), formula (12.351) can be deduced (see 
|42j ) from the Hamilton-Cayley theorem for symmetric matrices (cf e.g. [43]). 
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A general form of a quadrupole symmetric bitensor fielcfj], M 2j2 , can be found from 



fl2~35|) and ljAl7jl : 



2 2A 2A 

M 2 ^(a 2 ) = ^^mi 2A )(/3,C) £ (2^|2Air)r 2 ( ^ (2.37) 

A=0 n=A if =-2A 

(2\) 

with the six arbitrary scalar functions m„ . Since the fundamental tensors have the 
nonvanishing components only with even projections in the intrinsic system, the intrinsic 
components Mffl 2 v with odd |/i — v\ are equal to zero. Also, by virtue of (12.28al) and 
(12.2851) . the components differing in signs of both projections are equal to each other. In 
consequence, an arbitrary symmetric quadrupole bitensor has only six, and not fifteen, 
independent nonvanishing intrinsic components. It is convenient to work with their 
following combinations: 

M x = -MQ, - M$l_ x M = M 2 ( g 

M y = +M 2 ( g 1 - M^_! M, = V2M£ % (2.38) 

M z = -M { 2 \% + Mg_2 M 2 = M£ 2 % + Mg.2- 

These combinations are expressed in terms of some scalar functions m , m 1? m 2 , m' 2 , 
rri3 and m 4 (cf [6]) in the following way: 

M u = m — miP cos j u — m 2 (3 2 cos 2j u for u = x,y,z (2.39a) 
M = mo + m' 2 P 2 + m 3 p 3 cos 37 + (mi + m?,(3 2 )(3 cos 7 

+ (m 2 + m' 2 )(3 2 cos 27 + m 4 /3 4 cos 2 27 (2.396) 
Mi = — (mi + m.3/3 2 )/3sin7 + (m 2 + m 2 )/? 2 sin27 

— m 4) 5 4 sin 27 cos 27 (2.39c) 
M 2 = m + m' 2( 5 2 + m 3 /5 3 cos 37 — (mi + m 3 (3 2 )/3 cos 7 

- (m 2 + m' 2 )(3 2 cos 27 + m 4 /3 4 sin 2 27 (2.39(f) 

where the scalar functions appearing in ( I2.39al) . (12.3961) . (12.39cl) and (I2.39dl ) are related 
to the original scalar functions m?' in ( 12.371) by 



m 



mi 



m 2 



-j= (m< 4) /? 2 + 7t4 4) /? 3 cos37 + m 4 4) )^) (2.40a) 
mf(3 2 + 2mi 4) /5 3 cos 3 7 ) J (2.406) 

- mf/3 2 ) J (2.40c) 




m 2 = —rn K 2 3 m 3 = -m|, 4 ' m 4 = m^. (2.40c?) 



I A quadrupole nonsymmetric bitensor field is determined by specifying seven, and not six, arbitrary 
scalar functions. For such a field 
to the right-hand side of (T07j) . 



scalar functions. For such a field one more term, namely m\ (/?, £)(2{i2v\3K)x3K(ot2), should be added 
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Notice that because of the explicit dependence on (3 and 7 given in (12.3961) to (I2.39J) 
the six combination of the intrinsic components of M defined in fl2.38j) can be related 
to each other at some specific values of (3 and 7. For instance, at (3 — one has 
M x = M y = M z = M = M. Mi = 0, and at 7 = and arbitrary (3 it is M x = M y , 
M = M 2 and Mi = 0. 

3. The Bohr collective Hamiltonian 

The present-day notion of the Bohr Hamiltonian is not very precise. It encompasses 
a large class of Hamiltonians of which the original Bohr Hamiltonian [4] is only 
a very special case. Here, the general Bohr Hamiltonian means a generic second 
order differential Hermitian operator in the Hilbert space of functions of quadrupole 
coordinates oi 2ll , /•* = — 2, . . . , 2. Hamiltonians of a similar type but in the spaces 
of other collective coordinates, such as, for instance, the octupole deformations [32] or 
the pairing variables [44], can be called Bohr Hamiltonians as well, but they are not 
considered in this review. 

3.1. General form of the Hamiltonian 

The collective model arose from a classical description of nuclear collective phenomena 
(cf [H0]). This is why a classical Hamiltonian is often a starting point for its formulation. 
Now we follow this approach and construct the most general collective Hamiltonian 
using the quadrupole coordinates and making some natural assumptions. We assume 
the classical Hamiltonian to be: 

(i) a real function of coordinates a 2tl and velocities a 2fl 

(ii) invariant under orthogonal transformations of the coordinate system (the 0(3) 
scalar) 

(iii) an isotropic field (does not contain material tensors) 

(iv) a positive-definite quadratic form in the real and imaginary part of velocities a 2fJj 
and b 2lI . 

A general form of such a Hamiltonian reads 

H d (a 2 , d 2 ) = - ^B 2 fl /(a 2 )a 2u + V d (a 2 ) 

= - (a 2 ■ B 2 , 2 [a 2 ) ■ a 2 ) + V d (a 2 ). (3.1) 

The potential V c \(a 2 ) is a real and isotropic function of scalars do and xo- The inertial or 
mass bitensor B 2 2 is an isotropic symmetric field of the coordinates a 2fl . The bitensor 
matrix is such that the kinetic energy is positive. 

The quantum Hamiltonian corresponding to the classical one of (13.11) is obtained 
by the Podolsky-Pauli prescription [45l [46]: 

H coU (a 2 , ix 2 ) = — y~] 7r 2 M \f G(a 2 )G 2 ^ 2v (a 2 )i{ 2 + V co \\(a 2 ) (3.2) 

2a/G a 2 ~~ 
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where G = | det(i?2.2)| and G2,2 is the inverse matrix of the inertial matrix -82,2: 

52 B *V G *>*i = S v- ( 3 - 3 ) 

The Hamiltonian H co a (13.21) is Hermitian with the weight \fG. Hofmann [47j proved that 
the differential part of the Hamiltonian (13 .21) is unique provided it has the properties 
which will be specified below for the quantum Hamiltonian -ff qU ant, and fulfills the 
correspondence principle to the classical kinetic energy. However, the quantum kinetic 
energy is given only up to an additive scalar function (cf also [48]). This means that 
the quantum potential V co n in ( 13.21) is, in general, not equal to the classical potential V c \ 
of (|3.1I) and can contain quantum corrections. Usually, it is tacitly assumed that the 
potentials before and after the quantization are related by: V co \\(a 2 ) = V c \(a 2 ) + const, 
but this does not have any reasonable justification. The inertial bitensor £> 2 ,2 and its 
inverse G 2y2 are examples of symmetric matrices M 2fl ^ 2u (12.371 ). whose properties have 
been already discussed in section [2731 Thus, they are defined by specifying the six scalar 
functions. 

There are two ways to express H co \\ in the intrinsic coordinates. One way is to 
convert the derivatives djda 2il into derivatives with respect to the intrinsic coordinates 
using ( 127141) . (1272TD (and, if necessary, also (I27TTD . (127201) ( 127221 )1. The other way is to 
express the classical Hamiltonian H d (13.11) in terms of the intrinsic variables transforming 
the velocities using (I2.23P and (12.251) . One then has 

H c i = H c ^ TOt + -ff c i,vib (3.4) 

where 

= \ (409, iH + W, iH + W, lH) (3-5) 
-ffci.vib = 2 { B o( a o, 02)^0 + ZB^ao, a 2 )a a 2 + B 2 (a , a 2 )a 2 2 ) + V d (a Q , Xo) 

= X - ( J B w (A 7 )/? 2 + 2^ 7 (/?,7)/3/?7 + 5 77 (A7)/9 2 7 2 ) +Ki(AO- (3.6) 

The vibrational inertial functions Bpp, 5g 7 , -B 77 and B , Bi, B 2 are interrelated in the 
following way: 

Bppip, 7) = B {(3, 7) cos 2 7 + 2B t ((3, 7) sin 7 cos 7 + B 2 {f3, 7) sin 2 7 (3.7a) 

Bfrtf, 7) = (B 2 (P, 7) - B ((3, 7)) sin 7 cos 7 + B x {(3, 7 )(cos 2 7 - sin 2 7) (3.76) 

£ 77 (/3, 7) = B (J3, 7) sin 2 7 - 2B 1 (/3, 7) sin 7 cos 7 + B 2 (J3, 7) cos 2 7 (3.7c) 
and the moments of inertia read: 

4(00,02) = 4B u (a , a 2 )a 2 u = AB u (f3^)(3 2 sin 2 7^ ioxu = x,y,z (3.8) 

The functions B , B ll B 2 and B x , B yi B z are the corresponding combinations of the 
intrinsic components of inertial bitensor -B 2 ,2 (cf (12.380 ). Expressed in this way, the 
classical Hamiltonian (13.41) can be quantized with the Podolsky-Pauli prescription (cf 

EOT). 
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Both ways lead to the collective Hamiltonian expressed in terms of the intrinsic 
coordinates which consists of two parts: 

#coll = #rot + #vib (3.9) 

The rotational Hamiltonian H TOt is given by 

Hmt ~2 ^ WW) ( } 

u=x,y,z v-^' u 

and the vibrational Hamiltonian H vi ^ has the form 

H vih = h2 1 ( — xVg^— + Ij^Al 

2VGX \da G vlh da da 2 G vih da 2 

_ Ax^/- - ^-X^/G^-^-) + KoiiK, Xo) (3.11) 
oa Q G vib oa 2 oa 2 G vib o>a / 

k2 fl V^^+ „ 1 j- sin 3^-^ ^ 



2v / GV/^ 4 ^ G vib <9/? /3 2 sin3 7 <9 7 ' G v ib <9 7 

- Ttlk^T^^ ~ T^^T^^T^^) + ^(AC) (3-12) 
/3 4 d/3 G vib cry P sin 3 7 (9 7 G vib a/3 / 

with G vib = fio5 2 - B\ = BppByy - Bj 7 , G = B x B y B z G vih and X given by 

If the correspondence principle is not assumed i.e. if the considered nuclear 
collective system does not have a classical counterpart, one can construct a bit more 
general quantum Bohr Hamiltonian H quant possessing the following properties: 

(i) .ffquant is a second order differential operator in coordinates a 2 ^ possessing the finite 
lowest eigenvalue 

(ii) i/quant is a real operator (invariant under the time reversal) i.e. H qVLimt — -Equant 

(iii) Equant is a scalar operator with respect to the rotation group 0(3) 

(iv) -ffquant is Hermitian with a weight W(a , Xo) > 0. 

Under above assumptions iY qua nt can always be presented in the forrrj§|: 

Equant = -^77? ^2 ' A 2^2 ' M^2u + VquantC^O, Xo) (3.13) 

where A 2fl)2u is a symmetric positive-definite bitensor matrix which does not need to be 
related to W. Converted into the intrinsic coordinates with the help of (12.141) and (12.211) 
and pTB . (E2QD (ET22J), the Hamiltonian H quwt takes again the form with the 
rotational part ( 13.101) . The moments of inertia are expressed by the intrinsic components 
A x , A y and A z of A 2j2 (see (12. 38ft ) as follows: 

4a 2 4 

4(a ,a 2 ) = . - — r = -T-T-. — r/3 2 sin 2 7u for u = x, y, z. (3.14) 
A u (a ,a 2 ) A u ((3,y) 

§ When assumption (ii) is given up (reality of the Hamiltonian is not demanded) one more term can 
be added on the right-hand side of (|3. 13|1 . namely (1/W)(jr 2 ■ W / T 2 (a2)) + (T 2 (a2) ■ ^2), where T 2 is an 
arbitrary quadrupole field (cf [22]). 
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The weight W and the remaining intrinsic components of ^2,2, namely A , A\ and A 2 
enter the vibrational Hamiltonian in the following way: 

h 2 ( d d d d 

vlb 2WX \da ° da da 2 2 da 2 

d d d d \ 

+ —XWA,— + XWA,— 
oa oa 2 oa 2 oa J 

+ V r qU ant(crO,Xo) 

+V qmrt (/3,C) (3.15) 

where the relations between App, Ap 1 , A 77 and A , Ax, A 2 are identical with those 
between the vibrational inertial functions given in (13.7al) to (I3.7d l. 

The Hamiltonian (13 .21) obtained from its classical counterpart (13 . 1 1) by the Podolsky- 
Pauli quantization procedure is a special version of (13.130 in which W = vG and 
^2,2 = G 2 ,2- In this case W and ^2,2 are related to each other and both come from the 
classical inertial bitensor -82,2- Below, to the end of section [3J the Hamiltonians H co \\ 
and i/quant will not be distinguished from one another. Both of them will be denoted 
simply by H and their potential parts by V. 



3.2. Collective multipole operators 

Nuclear collective states are most often investigated experimentally by means of gamma 
spectroscopy. The collective model should therefore make it possible to calculate 
measured quantities such as the electromagnetic transition probabilities as well as 
the spectroscopic and intrinsic moments. To this end the electromagnetic multipole 
operators in the collective space should be constructed. It is assumed that the 
electric multipole operators are isotropic tensor fields dependent only on the collective 
coordinates a 2 ^. Since the fields depending on a 2 all have the positive parity, only the 
electric multipole operators with even multipolarities can be constructed. The electric 
multipole operator with an even multipolarity A is proportional to the 'dimensionless' 
A-pole moment of the charge distribution, Q^ ha,rge \ and can be written in the following 



fornijj: 



M(E0) =-=ZeI%qW(a ,Xo) (3.16a) 

V47T 



for A = and 



M(EX;n) = ^ 2 -^ZeR*Q^\a 2 ) 



|| The electric multipole operators are not the differential operators. They are denoted with the hat 
just to keep the notation for all electromagnetic operators uniform. 
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2ATT a/2+1 

— — ZeRl S\ ql X \o- ,Xo)Tk\n{a2,o- 2 ) (3.166) 

l07T * — ' 

fc=l 

for A = 2, 4, .... In these formulae the functions gj^ are scalar functions and 
R = tqA 1 / 3 . The fundamental tensors r fcA appearing in (13.1661) are given by (12.30d) . 
The intrinsic components of the moments of charge distribution are given by 

A/2+1 

Qit rgC)(in) (/5,7) = £ ^(/5,C)rg(/5,7). (3.17) 
fe=i 

According to (I2.31l) - (12.33l) the non-zero intrinsic components of the quadrupole and 
hexadecapole moments read: 

Q?o harge)(in) (/5,7) =g! 2) /3cos 7 + g( 2) /5 2 cos2 7 (3.18a) 

Q^ Tge)im \^l) = -^(<?i 2) /?sin 7 - gf/? 2 sin2 7 ). (3.186) 

g (charge)(iB) (/3)7) = 1 (? (4)^2 (5 ^2 ? + 1} + 2q^ f3^ COS 7 ( 7 COS 2 7 — 4) 



70 

+g$ 4) /3 4 (5cos 2 2 7 + 1)) (3.19a) 



g (charge)(iB) (/3)7) = I ^( 2g ( V sin 7 cos 7 - gf /3 3 sin 7 - 2gf /3 4 sin 2 7 cos 2 7 ) (3.196) 

g (cha rg e)(in) ( ^ 7) = 1 (g (4 )/?2 ^2 7 _ 2g (4)^3 ^2 ? ^ ? + ? (4)^4 ^2 ^ ( g ^ 

(2) (2) 

For the E2 moment it is often assumed that, in approximation, q[ = const and q 2 = 0. 

Similarly, only the magnetic multipole operators with the odd multipolarities can 
be constructed. It is assumed that they all are the odd-rank tensor fields dependent 
linearly on the angular momentum operator ( 12.21) . They can be written in the form 

A+l 



/2A 4- 1 Z 1 
M(MA; fi) = J ——-/IkR*- 1 £ [T K (a 2 ) x -L^ (3.20) 

71 K=A-1 

for A = 1, 3, . . ., where, according to (12.351) , the gyromagnetic tensor fields T K for k > 2 
are given by 

k K 

r KM (a 2 ) = E5 , I K) (o-o,Xo)^(«2,o"2,X3) (3.21) 
fc=i 

and T = g^\ao, Xo)- The functions gj^'s are called the scalar gyromagnetic functions. 
It follows from (I3.20f) that the intrinsic Cartesian components of the Ml operator read 

M U (M1) = y^^(/3, 7 )L„ (3.22a) 

9Ml) = 9o 0) ~ \J\g? ] ^os lu - y|^ 2) /3 2 cos2 % (3.226) 

for u = x, y, z (cf [7] for the j3- and 7 -dependence of g u ). It is often assumed that 
Qq^ = 1 and g± = g^ = (cf e.g. equation (1.37) in [49]). The coefficients in front 



Bohr collective Hamiltonian 



18 



of the right-hand sides of (13. 16 ah . (13.1651) and (13.201) are chosen to ensure the proper 
physical dimensions but their specific form is, to a large extent, a matter of convention. 

3.3. Collective wave functions 

In the previous sections the kinematics and dynamics of the collective model have been 
formulated. The model is determined up to a number of scalar functions defining the 
collective kinetic energy as well as the potential and the electromagnetic multipole 
operators. The question arises whether the collective wave functions can also be 
expressed through a number of scalar functions. The answer is 'yes'. However, the 
search for such scalar functions can turn out to be impractical. 

The collective wave functions are common eigenfunctions of the three operators, 
H, L 2 = (Li • Li) and L w , namely 

HVumM = Eu^umM (3.23a) 
L 2 V lIM (a 2 ) =h 2 I(I + l)V iIM (a 2 ) (3.236) 
L 10 ^ iIM (a 2 ) = hMV iIM (a 2 ) (3.23c) 

where i simply numbers (numerical) solutions or stands for a set of three additional 
quantum numbers (in the case of analytical solutions). Since equations ( 13.2361 ) and 
( 13.23 d) are automatically fulfilled by the wave functions of the form 

^umM = E Vkiijo'o, XoWimM (3.24) 

k 

where r^i are fundamental tensors ( 12.351) and vuu are arbitrary scalar functions, it only 
remains to determine the functions Vku using (13.23al) . In the intrinsic coordinates the 
wave function (13.241) has the form 

K=-I 
K k 

= E E MP> OtkiK^D^ip, t?, V)- (3.25) 

K k 

The exact analytical solutions of (13.23al ) are known for several special forms of 
the Hamiltonian (13.91) . In all of these cases the full inertial bitensor is replaced 
by only one constant mass parameter which (in view of fj2.39aD to (|2.39oT) ) leads to 



B x = By = B z = Bpp = i? 77 = B = const, B^ = (kinetic energy used originally by 
Bohr |4J). A comprehensive review of the exact and approximate solutions for different 
potentials is given in [50]. Here we discuss only a few cases with analytical solutions. 

We start with an obvious remark that if the potential has the form of V((3, () = 
Vp(P)+V y (C) I P 2 , the variable (3 can be separated from the remaining angular coordinates 
in the eigenvalue equations (13.23al) . (13.2361) . ( 13.23d) . In consequence, the scalar factors 
in the wave function (13 .251) factorize as follows 

<f> knAvI (P,0 = ^(P)wtm ■ (3-26) 
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Here the index i is replaced with the set of quantum numbers (nAu) with n labelling 
the energy levels E n &, A being here the separation constant, and v — an additional 
quantum number introduced by Arima (cf [51]) to distinguish different solutions with a 
given A. The variable (3 plays the role of the radial coordinate in the five-dimensional 
space and stands for the corresponding radial wave function. 

The angular part of the wave function (13.251) 

y&( 7> = «i?(0M7)^, VO (3.27) 



A' 



fulfills the equation 

sin 37 #7 07 at 4 sin 7„ at y 

= AY%( r ,,<p,0,iJ;). (3.28) 

The most important and most extensively studied is the case with Vy(Q =0, in which 
the separation constant A is equal to A(A + 3), where A, called the seniority [52], is a 
natural number. In this case a set of kj ordinary second order differential equations for 
functions (k — 1, . . . ,kj) with a given A and I = 2L + 3s (s = 0, 1 for even, odd 
J, respectively) was derived in [53l l40l [42]: 

9 _d (1 _ ^ ^fcWC) _Q {L + k + 3s)C ^k2L + 3s(0 

+ (A — 3s — L — k)(\ + 3s + L + k + 3)w£ L+3s (() + 4(k + l)(k + 2)w& 22L+3s (() 
+6{L _ k + ^ M-^L+ssiO + u{k + 1) dt< 1 a^a.(0 = Q (3 2g) 

(z/ labels the possible different solutions). The system of equations A3.29H can be solved 
by means of the method of polynomials. Bes solved this system for / < 6, A < 9 in [53]. 
Polynomials w£j for all values of the quantum numbers were found in [40], [54] using the 
theory of harmonic homogeneous polynomials. The functions Y^(j, ip, ip) form a set 
of spherical harmonics in the five-dimensional coordinate space (cf [55]). Another set of 
the five-dimensional spherical harmonics as functions of the biharmonic coordinates in 
the laboratory frame was constructed in [56]. Moreover, note that the operator in the 
l.h.s. of (13.28p is proportional to the Casimir operator of the SO (5) group and algebraic 
methods can be applied to calculate the SO (5) Clebsch-Gordan coefficients and reduced 
matrix elements of tensor operators, see [29} [30, E7[ [58] . Let us add that the solutions 
of (13.281) for some simple nonzero Vy(C) potentials are discussed e.g. in [59~l l60[ l3l]. 
The radial equation reads 

-j^h + £ + - u ' x,3) = °- (3 - 30) 

The elementary quantum mechanics (cf e.g. [61]) provides us with a few relevant solvable 
potentials Vp{f3) and we discuss some of them below. 

The Davidson modified oscillator potential Vg(/3) = \C{(3 2 + f3*J ' (3 2 ) has been 
applied for the first time in [62] to solve analytically the Wilets-Jean model [9]. Later 
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on the Davidson potential has been used by several authors (63J Ell EE]. The solution 
of the radial equation for this potential then reads: 

<<;<< » = J ,,. ' 1 \/^r i I v/^r r ) (3.31«) 





§ft(2n + «; + ^) (3.316) 
where Vn\x) are the Laguerre polynomials [41] and 

k(« + 3) = ^/5 c 4 q + A(A + 3). (3.32) 

The modified oscillator potential was also treated using algebraic methods in [3~T] . 

The Kratzer potential Vp(f3) = |C/?g q ((/9 e q//5) 2 — 2/5 eq //?) was introduced to the 
Bohr Hamiltonian in [66] ■ The Hamiltonian with this potential has a discrete spectrum 
only for negative energies. The radial functions and the energy levels of the bound states 
read: 

«£(/?) = N nX ^e- k ^L^\^-) (3.33a) 
=-^(f<)™ (3-33*) 



where k n \ = y2B\E n \\/h 2 , N n \ is a normalization constant, and k is given formally by 
(13321) . 

The Bohr model with the potential Vp(/3) in the form of an infinite square well has 
been discussed in [9], [67]. The corresponding radial wave functions and the energy levels 
are given by: 

u*(p) = N nX j^J x+3 /2(xn,\f3/f3w e \\) (3.34a) 
En> =^V^ (3-346) 

zn Pwell 

where N nX is a normalization factor, /3 we n is the well radius, J\+3/2(x) is the Bessel 
function [H] and x ni \ is its n-th zero. 

The modified oscillator potential for f3 eq =0 turns into the five-dimensional harmonic 
oscillator (the case of original Bohr Hamiltonian [4]), whose properties were studied most 
extensively. In this case the radial wave functions and the energy levels are still given by 
(13.31 aj) and (13.316|) . respectively, with k = X by virtue of (13.321) . For the angular part 



of the solutions one can take the functions 1/m(7> f, i>) discussed in the paragraph 
following (13.291) . The wave functions of the five-dimensional harmonic oscillator, but in 
a form which does not make their tensor structure explicit, were found also in [35] and 
[68j. Obviously, the eigenvalue problem for the harmonic oscillator Hamiltonian can be 
solved immediately in the laboratory coordinates. The most convenient way to present 
the solution is to introduce the notion of the quadrupole phonon. The phonon creation 



Bohr collective Hamiltonian 



21 



and annihilation operators, (5^ and j3 2 ^ , respectively, are defined as follows: 
~ t 1 f{BC\ 1/4 fBCY 1/4 \ 

K = — | ( -fp) «2, t - 1 y-jp j 7r 2 „ I (3.35a) 
BC\ 1/4 u (BC\' l, \ u \ 

—) <-\- W ) <)■ (3-355) 
The Hamiltonian expressed in terms of the phonon operators has the well known form: 

H = ■ AO + (3.36) 

and its eigenstates have the structure 



\N[c]IM) 



% X . . . X $ 



N 



|0> (3.37) 



JM 

where the symbol [c] stands for a given coupling scheme and |0) is the ground state 
annihilated by /3 2 M : /^lO) = 0. In section [2731 it was shown how to establish independent 
coupling schemes. However, to have the states with a definite seniority A one should 
introduce the so called traceless creation operators [38l 142] . The methods of group 
theory have been applied to calculate matrix elements of the multipole moment operators 
between the basis states of the group chain U(5)dSO(5)dSO(3) [69|. Besides the exact 
solutions of the harmonic oscillator presented above, the approximate harmonic solutions 
of the Bohr Hamiltonian with the collective potentials possessing well pronounced 
minimum for j3 = /3 eq 7^ are known from the very beginning of the Bohr collective 
model (see e.g. [42] for a survey of them). 

Special cases of the Bohr Hamiltonian for which equation (13.23d) possesses exact 
solutions are important for various reasons, but it is clear that in order to solve (13.23d) for 
Hamiltonians outside this class one needs to use numerical methods. Equation (13.23d) 
can be solved by transformation into a matrix eigenvalue problem using an appropriate 
truncated basis in the collective Hilbert space or by direct numerical methods suitable 
for the second order partial differential equations. In both approaches one must take 
into account the specific properties of the functions belonging to the domain of the 
Bohr Hamiltonian. First of all, they must be invariant against the octahedral group 
transformations. Moreover, they must behave appropriately (it will be explained below) 
at the boundaries (including infinity) of the six sectors of the (/3, 7) plane discussed in 
section I2.2L 

The basis used to calculate the matrix of the Bohr Hamiltonian can be chosen 
as the set of eigenfunctions of one of the exactly solvable cases, e.g. of the harmonic 
oscillator which was employed by Dussel and Bes [701, Gneuss and Greiner [71], and Hess 
et al [72l Obviously, such bases automatically fulfill conditions mentioned in the 
previous paragraph. However, this does not need to be true for 'artificial' bases (used 
e.g. in [741 ESI EE [77], [78]) which can be constructed from an arbitrary complete set 
of functions which are square integrable with respect to the measure (12.11) . We briefly 
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discuss two examples of such a construction. The first one was proposed by Kumar [74] 
who, in order to ensure that the solutions have the correct tensor structure and satisfy 
the mentioned conditions, expands the scalar functions Vkn((3,C) °f A3. 25ft with given k 
and / in the following basis functions: 

V$ m \P, C) = e-^^-/ 3 o)V2 e - CTC (/5 3 C-/3o 3 Co)/2 (/5 2 _ /3 2)n (/3 3 C _ ( g gg) 

which depend on the four variational parameters: (3 , Co, &/3 and o^. These parameters 
obey the following auxiliary conditions: I = 1 and er<* = for flo = 0, and I = 2 for 
/3 7^ 0. In the second example another complete set of functions is adopted, this time 
in the space of functions of all five quadrupole variables /3, 7, <p, ip 



^?l(P,7,V,^) = e-^ 2 p n \ ™™l \D I MK {^) (3.39) 

I olll 1 1 i J I 

with n — 0, 1, . . ., m = n, n — 2, . . . , or 1 and —I < K < I. The modified, still non- 
orthogonal, basis after the appropriate symmetrization with respect to the octahedral 
group O of the transformations of the intrinsic frame (cf (13.43 aft . (13.4351) and (13 .43 eft ) 
has the following form 

= e-^n £ /M 2i(7) (DL^faM) + (-l)^_ 2/ (^^)) (3.40) 
1=0 

where functions JkjzA'Y) are some specific linear combinations of the trigonometric 
functions cosjry or sinjry with the possible values of p — m, m — 2, m — 4, m — 6 and 
I^(m) is the maximal value of Z for given I and m (see [77] and [78] for details). The 
functions (13.40ft are then numerically orthonormalized and the parameter a is chosen so 
as to make the basis optimal. 

Before we write down a system of partial differential equations for (13.25ft derived 
from (13.23aft and which can be solved by direct numerical methods [H [79l HO, EE], we 
discuss the consequences of the invariance of ^um with respect to the octahedral group 
transformations of the intrinsic system. This invariance leads to 

1 

% IM {P,l,ip^)= Yl ^k^lk)D{ IK { Vk ^ k ^ k ) (3.41) 

K=-I 

for k = 0, 1, 2, 3, where the intrinsic coordinates with the index k are given in table [3J 
It follows from (12T31) that 

1 

*SS(/?>7)= E ^U^lMA^Qk^k) (3.42) 

K'=-I 

For k — 0, this is, of course, the identity. Formula (13.421) implies for k = 1,2,3 the 
following properties of ^{™k 

*t&(Pn) = (-lY^lAPn) (3.43a) 
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*iS (P, 7) = (-±) K/2 ^?kP, -7) (3-436) 
*Sk&7) = E 4^(^7-2^/3)^,(0,^/2^/2). (3.43c) 

K'=-I 

From (I3.43a[) . (13.4351) it follows that = for odd J's and that the components with 
odd K all vanish. Moreover, it is seen from (I3.43al) . (13.4351) and (13.43 cl) that the value 
of the function tyfi±ic(f3, 7) for an arbitrary 7 (—it < 7 < 7r) can be determined by the 
values of \l/-^,(/5, 7) with different even, nonnegative values of K' and the corresponding 
7 from the sector < 7 < 7r/3 shown in figured! Thus, it is sufficient to find \l/^(/3, 7) 
with even K > in the sector 0<7<7r/3to determine the wave function for all values 
of 7. Finally, the wave function (13.251) can be presented in the form 

[I /"A ^(m) ,n \ 

% IM ((3, 7, <p, 0, VO = J2 Tl , ( D M2 k (^ tf) + (-l)^-a(^ V>)) (3-44) 

where [J/2] stands for the integer part of 1/2. If the wave functions ^u m (a2) are 
normalized to unity 

(i'I'M'\UM) = j % rM ,(a 2 )%iMMWdn = SnSnS^M (3.45) 
the normalization of the functions ^^(Z?, 7) is 

/•OO /-7T/3 I 7 / 2 ! 9 

6 / / E TX^'S^' 7)^09. 7) ^(/3 2 , /3 3 cos 3 7 )/3 4 sin 3 7 d/?d 7 

= (3.46) 
It is also convenient to introduce the intrinsic wave functions 



Mft^J [2/+ ;^ jM ) *S(ft7) (3.47) 

which can be referred to as the normalized intrinsic wave functions because by virtue of 
( 13461) and (13471) they obey: 

^2(^'i2k\^ii2k) = k'i (3.48) 

fc=0 

where the scalar product (<&i>i2k\®U2k) is meant as 

poo /*tt/3 

(Wl^i/a*) = 6 / / $*, m (/5,7)^/2 fc (/5,7)^(/3,cos3 7 )/3 4 sin37d/3d7 . (3.49) 
Jo Jo 

The wave function ( 13440 is expressed in terms of the normalized intrinsic wave functions 
as 

^ 7M (/3,7,<^,^) = 

2iTT [I/2] 

~^-T E "WA 7)Vl + 5feo (^Wv. V>) + (-l) J £^_ 2fe (<A 0, V)) (3-50) 

fc=0 
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A set of the differential equations for $u2k obtained from (13.23al) is the following: 



i7 vib (/3, 7) -Eu + ± (7(7 + 1) - (2A;) 2 ) U ' 



(2kf 



1 



,(/3, 7 ) I y ((3n)J 2I z ((3,j) 



U2k+2 



(3.51) 
(3.52) 



4 (/?, 7 ) 4 (A 7), 
where k = 0(1), . . . , 7/2(7/2 - 1/2) for even(odd) 7 and 

C /fe = V(l + 4i)(7 + 2k)(I + 2k - 1)(7 - 2k + 1)(7 -2k + 2) 

By virtue of the symmetry conditions (13.4361) and ( 13 .43 cj ) and of equations f ]3 . 5 If) the 
normalized intrinsic wave functions <&u2k(P, 7) fulfill the following boundary conditions 

(i) At 7 = 0: 



(ii) At 7 = tt/3: 



c?7 

$i/2 fc (/3,7r/3) 
^oo(/?,7) 



7=0 



$i/ 2 fc(/3,0) = forA^O. 

= for odd 7 
= 



(3.53) 

(3.54a) 
(3.546) 



7=7r/3 



(7(7+l)-(2£0 2 )$, m (/3,7r/3) 

(/W3) + C /fe $ i/2fc _ 2 (/W3)) 
(7(7 + 1) - (2A;) 2 - 8) (W^, 7 )/a 7 )| T=7r/3 

' CWi (<9$ 

U2k+2 (A 7)/^) I 



> for even 7. 



l7=7r/3 



+C /fe (^ m _ 2 (/5, 7 )/9 7 )| 7=V3 



At /3 = the boundary conditions are 

d$ iI2k ([3,j) 

13=0 



d/3 



for 7 = 



and 



(3.54c) 



*i/2fc(/3 = 0,7) = for 7^0. (3.55) 



The bound-state wave function ^u2k(P,l) should vanish at infinity (J3 — > oo) faster 
than 1//3 2 to ensure its norm be finite. Equations (I3.5ip to (I3.55P are the basic ones for 
direct numerical calculations of <&u2k- Note, however, that such calculations are quite 
difficult due to the complexity of the system of equations (I3.5ip (in addition, the number 
of equations grows as 7/2) and of the boundary conditions. 

We end this subsection with the formulae for the reduced matrix elements of the 
collective electromagnetic multipole operators. From ( 13.1651) . (13.171 ) and fj3.50h the 
reduced matrix element of EA is 



<z'7'||M(EA)||z7) 



(2A+ 1)(27+ 1) ZeRl 



16tt 2 

x + M(l + S kl0 )(I2kX K \I'2k')(^ r2kl \Q^^ a) 

k,k'.K 



(3.56) 
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Hence, the E2 reduced matrix element reads 



167T 

fc>0 



■U2k) 



+ V 7 ^ 7 ^; ((/2fc22|r2fc + 2)($^ 7 , 2fc+2 |gf 2 hargc)(iri) |^ /2& ) 

+ (I2k + 22 - 2|/'2A;)(<|., r2fe |gf 2 harge){in) |<|. im+2 )) } . (3.57) 

The general formula for the MA reduced matrix elements is more complicated. We give 
only the formula for the Ml reduced matrix element: 



(j'I>\\M(Ml)\\U) = J m +1) ^ifJ2 {2k(I2klO\I'2k)(^ r2k \g z \^ iI2k ) 



An Ah 

k>0 



' ('^k)(I-2k + l) (I2k _ n \ 

I ~ 2k)i { + 2k ^(I2k + Ul\I'2k + 2)($,, P2M \ 9 -±-^\$„ 2t ) 



(I + 2k + 2){I 2k l) (m + n _ 1|m)(t ,_, / , 2t| « L _J ! , |$im+2) | ^ (3 M) 
5.^. Collective sum rules 

Quantities measurable experimentally are expressed, generally speaking, through matrix 
elements of the observables Ta^o^j 712)5 namely through the quantities 

(i'I'M'\f Xl ,\iIM) = J (^t'rM'Mf x ^a 2 ,w 2 )% IM (a 2 )^Wdn. (3.59) 

By comparing matrix elements calculated in the frame of a theoretical model with 
corresponding experimental values one can judge the quality of the model. Some 
restricted tests of the correctness of the theory are provided by the so called sum rules. 
The notion of sum rules is comprehensive and not too rigorously defined. The ones which 
we discuss below are called the collective sum rules, because they are restricted to the 
space of collective states (cf [82]). First we derive the so called non-energy weighted 
sum rules [83]. They were proposed by Kumar [84] to determine the intrinsic electric 
quadrupole moments of the nuclear states from the experimental data. Further on, the 
energy weighted sum rules [85] will be discussed as well. 
From the closure equation in the collective space: 

^2\iIM)(iIM\ = l (3.60) 

iJ.M 
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and the Wigner-Eckart theorem [34] the mean value of a tensor field T x (a 2 ) squared can 
be written as a sum of products of the reduced matrix elements of T A 

(™|(T A -T A )ii/M) = ^^^i(i/||T A ||i'/')| 2 . (3.61) 

i'l' 

Similarly, for the expectation value of the scalar of the third order in T A we have 

( ~ 1)A (iIM\([T x x T\]\ ■ T x )\iIM) = (UM\[[T X x T A ] A x T x ) \iIM) 
V2A + 1 

" N X X (^ll^iiz!/!)^!/!!!^!!^)^^!!^!^/) . (3.62) 

Similar formulae for higher powers of T A also be derived. However, for the fourth and 
higher order scalars in T A it becomes essential that the field T x should not depend on 
the momentum 7r 2 and its components should commute with each other (cf [83]). If 
it is possible to measure in experiment sufficiently large number of the reduced matrix 
elements of T A , equations f)3.6ip . (13.621) make it possible to determine the quantities 
on their left hand sides. Some of these quantities are otherwise difficult to access 
experimentally. When T 2 = M(E2) formula (13.611) reads 

(UM\ (M(E2) • M(E2)) \UM) = ^ B(E2; il -> i'l') (3.63) 

i'l' 

which means that the mean value in the state \iIM) of the intrinsic electric quadrupole 
moment squared is equal to the sum of the reduced probabilities of the E2 transitions 
to all accessible states \i'I'M'). The sum rules for the E2 operator have been applied to 
the experimental determination of the intrinsic quadrupole moments in |82j. 

The so called energy weighted collective sum rules are derived from the following 
double commutation relation: 

= K 2 A 2lh2v (a 2 ) (3.64) 



«2 



Hi 



H, a 2u 



where H is the collective Hamiltonian (13 . 1 3[) . Taking the mean values of both sides of 
(I3.64p in the state \iIM) and making use of the Wigner-Eckart theorem one obtains 

h 2 (iI\\A 2X (a 2 )\\iI) ^ /rn f 2 2 2A 



2 V4A + 1 



^(Ei'i'-E^i-iy-n ; ; z f ; i(z/iia 2 ||z'/')i 2 ( 3 . 6 o) 



where 

2u(oc 2 ). (3.66) 



For A = the energy weighted sum rule (13.651) takes the form 

h 2 

g v lv --^ ■ ~-y ■ ■ -- u ■ , 2 J _|_ _ 

i',I> 

The energy weighted sum rules can therefore be used to estimate the average values of 
the inertial functions using the excitation energies and the reduced matrix elements of 
a 2 extracted from the experimental data [85]. 



-(iIM\{A x + A y + A z + A + A 2 )\iIM) = ^(£^ - ^/)|(2/||a 2 ||zT)| 2 .(3.67) 
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4. Microscopic theories of the Bohr Hamiltonian 

The present section focuses on the derivation of the Bohr collective Hamiltonian from 
the microscopic theory based on an effective nucleon-nucleon interaction. There are two 
main methods which can be used to achieve this aim. The first one is the adiabatic 
approximation of the Time Dependent Hartree-Fock-Bogolyubov theory (ATDHFB) 
and the second one is the Generator Coordinate Method with the Gaussian Overlap 
Approximation (GCM+GOA). Before discussing the details concerning the application 
to the Bohr Hamiltonian case, we briefly recall basic assumptions underlying these 
methods and their general results. 

In the following it is assumed that the nuclear many-body Hamiltonian is given by 

H micr = ^ Kpud+dy + - ^2 V ^apd+d+dpd a (4.1) 

with the two-body interaction term V that can depend on the nuclear density, while 
d+ and d v are the nucleon creation and annihilation operators. The methods presented 
below can be directly applied to most of the effective non-relativistic nucleon-nucleon 
interactions. Some aspects of the Relativistic Mean Field (RMF) theory need, however, 
some additional discussion. 

Within the mean field approach the ground state of a nucleus is described by a 
generalized product state (in other words a BCS-type state, with the Slater determinant 
as a special case). To describe collective phenomena one needs a family of such states 
which depend on properly chosen collective variables such as quadrupole variables 
discussed in the previous sections. In the next step the collective Hamiltonian is 
constructed as an operator in the Hilbert space of functions of the collective variables. 



4.1. The ATDHFB method 



The ATDHFB method is conveniently formulated by using the notion of a generalized 
density matrix which for a given BCS-type state \^) and a given set of creation and 
annihilation operators cf+ , d n is defined as follows 



K 




where 



K 
K 

n 
n 



ii 

mn 

12 

mn 

21 

mn 

22 
'mn 



Wn"dm\V) 

(y\d n d m \ty 
(¥|cW+*l> 



(4.2) 



(4.3a) 
(4.36) 
(4.3c) 
(4.3d) 



The dimension of the 1Z matrix is twice the dimension of the one-particle space. A 
rigorous (but very abstract) definition of a space wherein the density operator acts can 
be found e.g. in [86]. In this section we will call it the double space. One should note 
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that there exist some slightly different sign conventions used to define 1Z. The 1Z matrix 
can be expressed in terms of the coefficients of the Bogolyubov transformation which 
connects the d+ , d n operators with the quasiparticle creation and annihilation operators 
corresponding to the state, see [871 [T5] . 

The Time Dependent HFB equation of motion reads 

ihK=[W(K),K] (4.4) 

with W being the self-consistent Hamiltonian (in the double space) induced by the 
density matrix 1Z 

™, ( K ~ XI \ ( F A \ 

w ={ o -k* + xi ) + {-A* -r* ) (45) 

in which A is the chemical potential and 

r^v = 22 V w'w'Pv'n' + rearrangement terms (4.6) 

= - V^v^viK^vi . (4.7) 

The rearrangement terms in 04.61) appear if i7 micr depends on the density. Typically 
these terms are considered only for the density dependent interaction in the particle- 
hole channel. 

In the adiabatic approximation the 1Z matrix is expressed as 1Z = 
exp(ix(£))7?-o(t) exp(— = TZq + TZx + TZ^. + ■ ■ ., where x is a 'small' operator, cf 
[T6] . The W operator can then be written as W = W + Wi + W 2 + . . ., with W fe 
induced by the TZk terms. Note that for k > only the second term in (I4.5P is used. 
By inserting the expansions of 7Z and W into (I4.4H one obtains a sequence of equations 
from which only the first two are usually considered 

ihK =\W ,K 1 ] + \W 1 ,K ] (4.8) 

inn, = [Wo, n 2 ] + [w x , + [w 2 , n ] + [w , n ] . (4.9) 

From the second equation one can infer that [Wo,7?.o] must be of the second order in 
the smallness parameter. An important feature of the r.h.s of the first equation (14.81) is 
its linear dependence on the TZi matrix. 

In the case when 7Z depends on time through several collective variables a. = 
(ai,... ,a n ) equation ( 14.81) leads to 
dlZ 

ihd k — - = [W ,7tf] + [Wi(7^),72o], k = l,...,n (4.10) 
where, since TZi depends linearly on 7Z , 

IZ^J^Ki- ( 4 - n ) 
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The collective classical Hamiltonian is defined as the expectation value of the H micY 
Hamiltonian in the |^(a)) state and can be written as a sum of the kinetic and potential 
parts 

H cl = T cl + V cl (4.12) 

where 

V cl = (V (<x)\H micr \V ( a )) (4.13) 
with |\&o) being the BCS-type state corresponding to the IZq matrix and 

k,j 

The inertial functions (also called mass parameters) B k j read 

B kj (ct) = ^Tr 2d (ni [||^o]) (4.15) 

where Tr 2( 2 denotes the trace in the double space. It can be verified that B k j = Bj k . The 
positive definiteness of the matrix B is not shown in the most general case, however it 
can be proved for several important instances, e.g. for cranking inertial functions, see 
[87] . The classical expression for H d needs to be quantized. Usually this is done by 
taking as the quantum mechanical kinetic energy the Laplace-Beltrami operator (times 
— h 2 ) in the space of the variables a. with B as the metric tensor and as the quantum 
potential energy the operator of multiplication by the function V c \ (this is the so called 
Podolsky-Pauli method of quantization, see also section l3TTi ). 

In the following it is assumed that (14.81) yields the unique solution for TZi in terms 
of 7^o • For a discussion of the existence and uniqueness of such a solution see [87]. 
However, one should remember that in practice the calculation of 1Z\ can be quite a 
difficult task. In the case of vibrational inertial functions only one collective variable 
(for axial symmetry systems) was considered, see [87] and [88] (for the case without 
pairing). In [87] TZ\ was obtained by matrix inversion while the method proposed in [88j 
(see also [89]) was based on the fact that (14.81 ) can be transformed into (for simplicity 
we consider here a single variable a) 

[W + Wi - 2aP, TZ + Tlx] = (4.16) 

where 

v = ^[d a n ,n ] . (4.17) 

Equation (14. 1 61) can be solved by a constrained HFB calculation. In the case of a 
rotation (and in the limit of vanishing angular velocity) it leads to the Thouless-Valatin 
expression for the moment of inertia [90]. For this reason, the term [Wi, TZo] is sometimes 
called the Thouless-Valatin (TV) correction. 
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The trace in (14.151) is most easily calculated in a quasiparticle basis (in the double 
space). It can be verified that both dj-TZo '■— dTZo/dak and TZ{ have the antidiagonal 
form 




(d k TZ ) qp = (4.18) 



z 



= n ( 4 - 19 ) 



so that 



Tr 2d {n{ [d k TZ , n }) = Tr^tf - z*f k ) . (4.20) 

Neglecting the TV term brings in a substantial simplification, because then 

i^ttfc/fc,^ = (E^ + E v )z ktllv (4.21) 

where are quasiparticle energies, and to calculate the inertial functions in this 
approximation one needs only d k 7Zo and E^. 

S * = T^ (^T^) • (422) 

Inertial functions calculated without considering the Thouless-Valatin term are 
sometimes called the cranking inertial functions. 

Matrix elements f k>fll/ can be expressed in terms of the derivatives of the density 
matrix p and the pairing tensor k in the canonical basis 

fk {d k p)^{u^v v + u M u„) + (dk^^u^Uv - VpV v ). (4.23) 

In this formula s v is the phase factor with the properties = 1, s M = — Sp_, while p 
and p denote canonically conjugated states. Another, perhaps more familiar, expression 
with the explicit dependence on the derivative of the BCS state is 

fk,^ = (*o|a„tv|d*tf ) (4-24) 

where a M are the quasiparticle annihilation operators. If the condition [Wo,7£o] ~ is 
used one gets yet another, sometimes quite useful, form 

fk,iw = ~ F 1 F [su(dM^(u^v u + VpUv) + (d k A)^(u^u u - vflv)] (4.25) 

where h = K — XI + r, see (14.51) . The references cited in section 14.31 contain many 
examples of application of formulae (14.221) to ( 14.251 ) . 
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4.2. The GCM+GOA method 

This method is based on the variational principle. Its simplest version consists of taking 
as test functions the integrals f dacf(ac)$(a) with the BCS-type states |$(o:)) and the 
weights /(a) which are to be determined. The assumption of the Gaussian overlaps 
(GOA) 

($(a")l$K)> = exp(- 9kj{<*){a'l - a' k ){a'> - a' j )/2), c* = {ct" + a')/2 (4.26) 

k,j 

enables one to transform the Hill- Wheeler- Griffin equation into the Schrodinger-like, 
second order differential equation for the weight functions /(«) 

H GC Mf(<*) = Ef(a) . (4.27) 

Below, only some of the final results for the GCM collective Hamiltonian Hqcm are 
presented. Details can be found in e.g. [2T] (see also [49, 23]). A specific application to 
the Bohr Hamiltonian problem was studied in [22]. The GCM Hamiltonian reads 

Hqcm = Tqcm + Vqcm (4.28) 

where 

^GCM = £ ' < 4 ' 29 ' 

Note that in the case of the quadrupole collective variables the Hamiltonian (14.281) has 
the same structure as the Hamiltonian (13.131) discussed earlier, with -y/det g = W and 
-B G cm = A. The inverse of the inertial tensor is given by 

(Bohu)" 3 = ^( Re/ ^ - Re/m)« (4.30) 

where 

hn,mn = D a ^D a ^h(a" , a') | a "=a'=a (4.31a) 

hl2,mn = D , l i n D a ; n h{0L" , a!) \ a "=a!=a (4.316) 

with 

h{d',a!) = ($(a")l^mic r |$(a / ))/($(«' / )l$(« , )) ■ (4-32) 

D denotes here a covariant derivative and the indices on the right hand side of (14.301 ) 
are raised using the inverse of the metric tensor g^j. The tensor g k j can be written as 

fly(a) = (d afc $(a)|0 aj $(a)) . (4.33) 

The potential energy can be transformed into the form 

^gcm = V cl (a) + V ZPE (a) (4.34) 



where 



V^a) = ($(«)|# mi c r |$(a)> (4.35) 
*W«) = -% Y^^Gcydki ~ \ ^2g kj D ak D a .V c i(a) . (4.36) 



k,j k,j 
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The first term V c \(a) in (14.341) is essentially the same as in the ATDHFB method. The 
second one, Vzpe(ck), specific for the GCM method, is called the zero point energy 
correction. One of the advantages of the GCM over the ATDHFB method is that it 
avoids ambiguities of the quantization procedure. 

Let us now add two comments. Firstly, the derivation of (14.281) is based on the 
assumption that the space with the metric tensor g is essentially flat, which means that 
one can introduce such variables that g is constant. In some references, e.g. [91], it 
is assumed from the very beginning that collective variables do already have such a 
property. Secondly, the overlaps (14.261) can in principle contain imaginary terms linear 
in ex." — ex' , which would lead to extra first order derivative terms in (14.281) . see [21, 92J. 

It seems that the GCM formulae for mass parameters, which in the case of the Bohr 
Hamiltonian were discussed for the first time in [22j, have never been used in practical 
calculations. There are two possible reasons for this. The first one is their technical 
complexity as compared to the cranking inertial functions of the ATDHFB method. 
The second one is the known deficiency of the simplest version of the GCM, which can 
be easily seen in the case of the translational motion, for which the GCM fails to give 
the correct mass. This deficiency can be amended by considering a richer set of test 
functions with the so called conjugate parameters. The even more sophisticated version 
of the GCM presented in [93], EE] leads to the inertial functions which are the same as 
those in ATDHFB and to the potential energy containing the ZPE terms, see also [94J. 

4-3. Application to the Bohr Hamiltonian case 

Studies on the microscopic foundation of the Bohr Hamiltonian have a long history. 
In fact they were an important stimulus for the development of the ATDHFB method 
(see [6] and [15], where the general microscopic formulae for the inertial functions and 
potential energy entering the Bohr Hamiltonian were obtained). The first attempt to 
apply these results to the real nuclei (in the W-Pt region) and to compare them with 
experimental data was performed by Kumar and Baranger in [95], see also [74]. They 
used a spherical mean field potential plus a separable quadrupole-quadrupole interaction 
and the seniority (constant G) pairing. In the following years their quadrupole+pairing 
model has been applied to several regions of nuclei, to cite only the most recent papers 
[961 E7]. 

The second attempt, based on a picture of nucleons moving in a deformed single- 
particle Nilsson potential was presented in [98], [99]. The collective potential energy 
was calculated by means of macroscopic-microscopic method using the Strutinsky shell 
correction method. More recently this model has been extended by an approximate 
inclusion of the effects of pairing vibrations [77|, 11001 |82j. 

There are several technical difficulties, including a substantial computational effort, 
if one wishes to apply the theories presented in previous subsections to the case of a 
nuclear mean field obtained self-consistently from realistic effective interactions, e.g. of 
the Gogny or Skyrme type. The first results with the Gogny force (GCM+GOA) were 
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announced in |101| . but more extensive calculations were done only in the 90's |102l [78] . 
see also more recent |1 031 1104] . The Skyrme interaction was used in the framework of 
the ATDHFB in [76], and next, in a more extended way in [1051 11061 1107j . Let us also 
mention the papers |108U109j . where the Skyrme force (with GCM+GOA) was used with 
an interesting, although somewhat simplified treatment of the Bohr Hamiltonian. In 
addition, recent years brought applications of the RMF to describe quadrupole collective 
states |110l 1111] . The range of nuclei studied in the cited papers is quite wide: from the 
krypton isotopes to transactinides. 

An important advantage of the microscopic approaches is a clear physical 
interpretation of the quadrupole variables, which in most cases are simply proportional 
to the components of the quadrupole tensor of the mass distribution, that is in the 
laboratory frame 

a 2 ,~($l< aSS) | $ > (4-37) 

with 

^r S) = X>^(^) (4-38) 
i=l 

where the sum runs over all nucleons and A is the mass number. A slightly different 
model is the one studied in [98II99], where the collective variables describe an elliptic 
deformation of the single particle potential. In the intrinsic frame, instead of Q^^, 
one often uses the Q , Q 2 operators: 

Q = £(3z? - T?) = V^F/5QSr ) (4.39a) 

i=l 
A 

Q2 = Y. - Vi) ■ (4-396) 

i=l 

The intrinsic coordinates (ao,a 2 ) and 7) (see (12.5al) . (12.561 ). (12.61) ) are given by the 
mean values of the Qo, Q 2 operators 

a = /3cos7 = cq = c($\Q \$) (4.40a) 

a 2 = /3sin7 = cg 2 = c($|Q 2 |$) (4.406) 

with the coefficient c conventionally chosen as 

c = (4.41) 

to get some correspondence with phenomenological models (see e.g. [42]). For r 2 the 
liquid drop model estimate ^{tqA 1 ^) 2 with r = 1.2 fm, is usually adopted. 

All microscopic theories start with BCS-type states constructed in the intrinsic 
frame which are subsequently rotated to a desired orientation. In the self-consistent 
theory a set of states parametrized by the quadrupole variables go, <?2 (or equivalently 
P, 7) is obtained from the HFB calculations with constraints on the mean values of Qo 
and Q 2 operators 

6{$\H micr - A0Q0 - A 2 g 2 |$) = (4.42) 
($|Qo|$> = g , ($|Q 2 |$> = g 2 . (4.43) 
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For the fixed values of go, <?2 the states |$) in the intrinsic frame should be invariant 
against the discrete symmetry group D 2 C O, which is generated by rotations by ix 
around three orthogonal axes. This property is obvious if one recalls that such rotations 
do not change the values of go and g 2 . An extensive analysis of the consequences 
of this and other symmetries for various quantities of the mean field theory can be 
found in |112l 1113] . Equation (I4.42j) is solved numerically by expanding the BCS-type 
states in a properly constructed basis. The basic ingredient of such a construction are 
the eigenstates of the three-dimensional harmonic oscillator in Cartesian |114l llllj or 
cylindrical coordinates |115j . 

Several properties of the kinetic part of the Bohr Hamiltonian in the intrinsic frame, 
e.g. the block structure of the inertia tensor, can be inferred from the general theory 
— see formulae ( 13.51) . ( 13.61 ). However, in view of the approximations adopted, it is 
instructive to prove these properties directly using the definitions of inertial functions 
given in section 14.11 We sketch the proof for the ATDHFB inertial functions in the 
Appendix B, for the GCM+GOA case cf [22]. 

In the intrinsic frame the rotational part of the collective space is identical with 
the SO (3) group. For the tangent space (in the differential geometry sense) of the 
SO (3) group it is more convenient to use a basis determined by the group structure 
itself, consisting of elements of the Lie algebra (in other words the angular momentum 
operators), rather than the usual coordinate basis (e.g. derivatives with respect to the 
Euler angles). The rotational part of the Hamiltonian (see (3.4)) is then expressed in 
terms of the collective angular momentum operators and the moments of inertia. The 
dependence of the BCS-type state | <&(/?, 7, <p, ip)) on the Euler angles is explicitly 
known (see Appendix B), and if the TV term is neglected, one arrives at the Inglis- 
Belyaev formula for the moment of inertia 



where is the angular momentum operator in the Fock space. Note that the sum in 
(I4.44p is over the whole one-particle space basis. 

The vibrational parameters are a bit more complicated, even if the TV correction 
term is neglected. To find them one can use (14. 23j) but the derivatives of p and k have to 
be in this case calculated numerically, see papers [89, 105] , where parameters obtained in 
this way were called the M p masses. More frequently, another approximation is adopted, 
which consists in using (14.251 ) but neglecting the derivative of A. The derivatives of IZq 
with respect to the Lagrange multipliers A, can be then calculated as 



(E„ + E v ) 



(4.44) 



(4.45) 



and next the derivatives of qi with respect to Xj 
dqi {»\Qi\v){v\Qj\tj) , 




(4.46) 
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As a result, one gets the widely used formula for the vibrational inertial functions 



The indices /x, v in (14.441) . (14.461) . (I4.48P run over all one-particle basis states, but in 
many cases, due to symmetry properties, their range can be limited to half of the basis, 
so one should take care when comparing formulae in various papers. The transition to 
the inertial functions for (3 and 7 variables, entering e.g. (13 . 1 2fl is straightforward. Note 
also that the symmetry properties of these functions, which can be inferred from the 
discussion of the symmetric quadrupole bitensors in section 12.31 provide useful tests for 
correctness of numerical calculations. 

Let us make one remark about the two kinds of nucleons. As long as one does not 
consider the proton-neutron pairing, the state |$) (see (14.421) ) is a product of the BCS- 
type states for protons and neutrons and consequently the ATDHFB mass parameters 
are simply the sum of the corresponding parameters for protons and neutrons. 

Once the mass parameters and the potential energy are determined, the eigenvalue 
problem of the Bohr Hamiltonian is solved numerically yielding the collective wave 
functions and energies, which can be compared directly with experimental data. The 
wave functions are then used to calculate matrix elements of collective electromagnetic 
transition operators, primarily of the E2 operator, which in the microscopic approach is 
equal to 



with the sum over protons only, and has the structure discussed in section 13.21 On 
the experimental side, modern techniques provide a large amount of data about the 
electromagnetic transitions. In several cases not only transition probabilities are 
measured but also relative signs of the matrix elements are determined, see e.g. [82]. The 
correct interpretation of rich experimental data is challenging and provides a stringent 
test of theoretical models. 

A very general conclusion stemming from numerous papers on microscopic 
calculations of the Bohr Hamiltonian is that it leads to a good qualitative agreement with 
experiment, correctly reproducing the main features of collective spectra. These results 
are remarkable, especially if one remembers that they are obtained without any free 
parameters fitted to collective quantities. At the quantitative level one finds, however, 
several substantial discrepancies between the theory and measurement. For example, 
quite often the theoretical spectrum is 'stretched' compared to the experimental one, 
although the two are similar in their general patterns. 

There are several possible reasons for such discrepancies. A natural question 
is what is the optimal nucleon-nucleon effective interaction in both particle-hole and 
particle-particle channels. As it is easily seen in the case of the Skyrme forces, various 




(4.47) 



(4.48) 




(4.49) 
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parametrizations can lead to rather different potential energies. On the other hand the 
inertial functions are quite sensitive to the strength of the pairing interaction. There 
are also some open problems in the collective part of the theory. Let us mention two 
of them — the question of the Thouless-Valatin corrections and of the extension of the 
collective space by including variables describing pairing degrees of freedom. The TV 
correction in the case of the moments of inertia can be estimated by performing cranked 
HFB calculations as described in [78]. A study of some test cases led to the conclusion 
that on average the TV term increases the moment of inertia by 30%. To account for 
this effect the constant factor 1.3 was used in [75] to scale the moments obtained with 
the Inglis-Belyaev formula. Much less is known about the TV correction in the case of 
vibrational parameters. There exist only old papers discussing the axial case [881 157], 
where the reported effects, depending on the nucleon-nucleon interaction, were of a size 
comparable to those for moments of inertia. In our opinion, in order to respect the 
symmetry conditions, when simulating the TV corrections in the simplified approach, 
all inertial functions, not only the moments of inertia, should be rescaled by a constant 
multiplicative factor. Such a common factor, equal to 1.2, was used e.g. in |105| . The 
question of extending the collective space is discussed in section 14.51 

4-4- Example of microscopic calculations 

The theory presented in the previous sections is now applied for twelve heavy even-even 
nuclei 178 Hf — 200 Hg. The considered chain of nuclei provides a good example of a 
transition from deformed axial ( 178 Hf) to almost spherical ( 200 Hg) nuclei. Theoretical 
results were obtained using the SLy4 version of the Skyrme interaction [116] and the 
constant G (seniority) pairing. The inertial functions and the potential energy entering 
the Bohr Hamiltonian were calculated using the ATDHFB theory (the formulae (14. 13ft . 
( 14.441) . ( 14.471) ). The general outline of calculations is similar to that of pTOBj, where one 
can also find more about the basis we used and other numerical details. Here, however, 
we do not use the scaling factor discussed in the end of the previous subsection. The 
pairing strength was determined by comparing the 'experimental' pairing gaps from the 
five point formula [117] for Hf and W isotopes with the minimal quasiparticle energies 
calculated for the deformation corresponding to the minimum of the potential. To get 
the inertial functions and collective potential energy, the HF+BCS calculations with 
linear constraints on qo and q2 were performed for 155 points forming a regular grid in 
the sextant < (3 < 0.7, < 7 < 60°, with the distance between points equal to 0.05 
and 6° in the /3 and 7 direction, respectively. 

A sample of the results is presented below. Figures EHU show the collective potential 
energy relative to that of the spherical shape for three nuclei: 178 Hf, 190 Os and 200 Hg. 
The inertial functions Bpp, Bp 1 , £> 77 , B xy ^ z are plotted in figures [5] and [6] for the 190 Os 
nucleus only. They exhibit quite strong dependence on the (3, 7 variables and differ 
significantly from the functions adopted in most of the phenomenological models, in 
which Bj3p—B 7 y—B Xi y iZ —B and Bp^ = (i.e. in which there is only one constant 
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Figure 2. Potential energy (relative to that of the spherical shape) for the 178 Hf 
nucleus calculated using the SLy4 Skyrme interaction plus the constant G pairing. 



A 




Figure 3. Potential energy for the 190 Os nucleus, see also caption to figureEJ 

mass parameter B). Next, the selected theoretical energy levels and the E2 transition 
probabilities are compared with experimental data taken from |118j . The levels are 
labelled with their spin I and the number k for a given spin as Ik- Figure [7] shows the 
levels 2i and 4 X belonging to the ground state band and the levels 2 2 and 2 , which are 
often treated as bandheads of the quasi 7 and (3 bands. The reduced probabilities of two 
transitions within the g.s. band, namely 2i — > Oi, 4i — > 2 X as well as of two inter-band 
transitions 2 2 — > 0\ and 2 2 — > 2i are plotted in figure El Of course, the results shown 
here and those not shown deserve a more detailed discussion which we postpone to a 
subsequent publication. However, one should admit that an overall agreement between 
the theory and the experiment seen in figures [7] and [8] is quite remarkable. It must be 
stressed once more that the sole input for the calculations is the effective nucleon-nucleon 
interaction without any additional parameters, effective charges, etc. 

4-5. Extension of the collective space 

We end this section with a brief discussion of ideas concerning the extension of the 
collective space to variables connected with the pairing degrees of freedom. Some vague 
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Figure 4. Potential energy for the 200 Hg nucleus, see also caption to figure [21 
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Figure 5. Vibrational inertial functions Bpp, Bp 7 , B 11 (in ft 2 /MeV) for the 190 Os 
nucleus. 
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Figure 6. Rotational inertial functions i?„, u = x,y,z (in fi. 2 /MeV) for the 190 Os 
nucleus. 



suggestions about such possibility can be found already in [6], HH]. In the simplest 
version of the collective approach to pairing, appropriate for the constant G interaction, 
the pairing gap A, which enters the formulae for the v^, coefficients, is treated as 
a variable and not a solution of the BCS equations. In this way a family of states is 
parametrized by A and the ATDHFB [44] or the GCM+GOA [22], [92] methods allow to 
obtain the collective pairing Hamiltonian. The case of the state dependent pairing needs 
a slightly different definition of a collective variable instead of A, see e.g. |119| . Usually 
one also considers a second variable, the gauge angle, which is used to implement the 
projection on a fixed number of particles. An important consequence of the collective 
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Figure 7. Lowest collective energy levels for 178 Hf — 200 Hg nuclei. Upper panels: 
2i and 4i levels in the g.s. band. Lower panels: 22 and O2 levels. Experimental data 
taken from |118| . 





Figure 8. Selected B{E2) transition probabilities (in W.u.) in the 178 Hf — 
200 Hg nuclei. Upper panels: transitions within the ground state band, lower panels: 
22 — * 0g. s . and 2 2 — * 2 g . s .. Note the different scales in various panels. Experimental 
data taken from |118| . 



treatment of the pairing is that one gets probability distributions in the A variable 
(coming from eigenfunctions of the collective pairing Hamiltonian) instead of a single 
value of the pairing gap. It is well known that values of the inertial functions entering the 
Bohr Hamiltonian strongly depend on the value of the pairing gap. Hence it is important 
to estimate the influence of the extension of the collective space on the spectrum and 
other properties of the Hamiltonian. 
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If both quadrupole and pairing variables are combined, the resulting extended 
collective space is nine dimensional (four pairing variables for two kinds of nucleons). 
In principle one can use the methods from sections 1431 or [4721 to construct the collective 
Hamiltonian also in this case [119J. However, a full analysis of such Hamiltonian would 
be very difficult so that one must resort to some approximations [77 1 1119j . Firstly, the 
terms which explicitly couple quadrupole and pairing variables are neglected. Secondly, 
the ground state of the pairing part of the Hamiltonian is determined to get the 
probability distribution in the A variable (for both protons and neutrons). In [771 HOOj 
the inertial functions of the Bohr Hamiltonian were calculated using the value of the 
pairing gap corresponding to the maximum of the distribution while in [119] the inertial 
functions obtained for various values of A were averaged with the weight given by this 
distribution. Results of the cited papers show that the collective treatment of pairing 
can bring significant changes in the energy spectrum, but more extensive studies (also 
on more consistent methods of determining the pairing strength) are still needed. 

5. Summary 

The Bohr collective model was originally invented to describe quadrupole oscillations 
of nuclear surface in the spherical nuclei. It has been subsequently developed and 
successfully applied to the description of the quadrupole collective excitations in various 
deformed nuclei. The Bohr Hamiltonian can be treated either at the phenomenological 
level or can be derived from a microscopic many-body theory. In the present review 
we have leaned towards the latter treatment of the Bohr model. Microscopic theories 
give, as a rule, quite complicated collective Hamiltonians which are determined by a 
number of functions of the dynamical variables rather than by a number of parameters 
as it used to be in the phenomenological approaches. Therefore, it is important to 
study the most general under some natural conditions, form of the Bohr Hamiltonian. 
The basic assumption of the Bohr collective model is that the dynamical variables form 
a real electric quadrupole tensor. The detailed analysis of the isotropic tensor fields, 
presented in this review, allows us to study the tensor structure of various ingredients 
of the general Bohr collective model. Its specific version is determined by several scalar 
functions. The collective Hamiltonian is defined by the six scalar inertial functions, 
the scalar weight and the (obviously scalar) potential — eight scalar functions of two 
scalar variables altogether. To define other collective observables we still need a number 
of additional scalar functions. The knowledge of the universal tensor structure of the 
collective wave functions allows us to express them also through scalar functions specific 
for a definite eigenstate of a given Hamiltonian. 

The scalar functions defining the model can be calculated in the framework of a 
microscopic theory. In practice it is the components of the inertial bitensor rather 
than the scalar inertial functions themselves that are calculated. There are the two 
fundamental methods to obtain the Bohr Hamiltonian from the many-body theory. 
One is the Adiabatic Time Dependent Hartree-Fock-Bogolyubov approach which, since 
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it treats the collective motion as the classical motion of a wave packet, gives the 
classical and not the quantum Hamiltonian. Quantization of the classical Hamiltonian 
obtained in this method can lead to an ambiguity in the collective potential. The 
other method is provided by the Gaussian Overlap Approximation of the Generator 
Coordinate Method. It directly gives the quantum second order differential collective 
Hamiltonian (without the Gaussian Overlap Approximation the Generator Coordinate 
Methods leads to integral equations). Both methods mentioned above can be applied 
to different microscopic models. A whole range of the nuclear many-body theories 
are available: the phenomenological mean field models with different single-particle 
potentials such as the Nilsson or the Woods-Saxon ones, the self-consistent approaches 
with the Skyrme or the Gogny effective nucleon-nucleon forces and the Relativistic Mean 
Field theory. The presented example of application of the reviewed approach to the 178 Hf 
— 200 Hg nuclei shows that the Bohr Hamiltonian derived from the effective interaction 
is a very useful tool for explaining nuclear collective properties. In addition, having 
the Bohr Hamiltonian obtained from the microscopic theory, one can try to examine 
another question which recently attracted attention, that is of the possible foundations of 
phenomenological models inspired partly by the Interacting Boson Model [671 [66]. This 
would require, however, a more extensive and detailed discussion which goes beyond the 
scope of the present review. 
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Appendix A. Spherical tensors and tensor fields 

In the nuclear collective models the operands are most often tensors under the 0(3) 
group of orthogonal transformations in the three-dimensional physical space. The 
models employ tensor fields such as potentials, inertial functions, etc., which are usually 
isotropic functions of the dynamical variables. Here, in order to fix the notions and 
notation, we collect the definitions and properties of quantities and operations appearing 
in the theory of spherical tensors. We also briefly present the structure of the isotropic 
tensor fields, which is not so widely known. 

Appendix A.l. Spherical tensors 

A set of 2 A + 1 covariant components with projections (magnetic numbers) fi = 
—A, —A + 1, . . . , A — 1, A, representing an object (a c-number, an operator, etc.) a\ in 
a three-dimensional coordinate system U, forms the S0(3) irreducible spherical tensor 
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of rank A if 



n\., DfafaW)^ (A.l) 



where a' x , are the components representing the same object a\ in the coordinate system 



U' rotated with respect to U by the Euler angles (p, ip } and D*, Jcp, ip) are the 



Wigner functions (rotation matrices) [34] [36]. The Wigner functions defined in [34] 
are complex conjugate of those of Bohr and Mottelson [36]. In nuclear physics the 
convention of Bohr and Mottelson is most widely used and we adopt it in this work. 
We do not quote here various useful properties of the Wigner functions as they can be 
found easily for instance in [34]. Only integer ranks A are considered here. 

In the special case of a rotation by the Euler angle with $ = 0, ip = (i.e. of a 
rotation about z-axis) formula (|A.ll) reduces to 



«v = E e ~ iw <w^ ( A - 2 ) 

p- 

A tensor whose only nonvanishing components is a\ (i.e. = for fi ^ 0) is invariant 
under such rotations, and can be called the axially symmetric tensor [z is the symmetry 
axis). In general, a SO (3) tensor transforming under rotations as in (jA.lft . does not 
have a definite transformation rule under the inversion of the coordinate system, which 
supplements the rotations to the 0(3) group. However, every SO (3) tensor can be 
decomposed into the two terms 

«» = «™+a™ (A.3) 

called the proper (polar or electric) and the pseudo- (axial or magnetic) tensor, 
respectively, with the following transformation properties: 

off* = (-DM? (A.4) 



and 



agr=(-l) w «2? (A.5) 



where the double prime denotes the components of the corresponding tensors in the 
coordinate system U" whose axes are reflected with respect to those of U. The sign 
factors on the right-hand sides of (I A. 41) and (IA.5I) are the parities. Further below, the 



superscripts (E) or (M) will be omitted when the parity is not relevant. In general, the 
components are complex numbers. However, in physical application most frequently 
real spherical tensors appear whose components are determined by 2A + 1 real numbers 
up to a single phase factor. More precisely, the complex components of such tensors 
satisfy the following relations 

ak«e«W(_i)-M OA _ M . (a.6) 

Relations (IA.6I) are preserved by the transformations f lA.ll) . In the following the common 
phase is set to zero (5(A) = 0) as in the case of the standard spherical harmonics Y Xfl . 
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Equation (IA.6I ) suggests that the tensor ot\ can be presented in the following form: 

i (-iy 

«ao = «ao OLx^ = "y|( a AM + ihn) oc\-n = -j= (a A/ , - ibxfj) (A. 7) 

for \x > 0. Thus, the tensor ot\ is indeed determined by the 2A + 1 real numbers: the 
real and imaginary parts, (/i = 0, 1, . . . , A) and (fi — 1, . . . , A), respectively, of 
the complex components for fi > 0. 

Coupling two real tensors, a a and Xk, of ranks A and k and the zero phases in the 
following way 

[a x x x*]lm = {X/ikv 1 LM) ax^Xnu, (A. 8) 

fJbV 

where (\nnv\LM) is the Clebsch-Gordan coefficient, one obtains a real tensor of rank L, 
where |A — k\ < L < A + k, called an irreducible tensor product. Its phase is determined 
by the relation 

[«a x Xk\lm = (-l) L - M - X - K [a x x x*]l-m. (A.9) 
The scalar product of the two real tensors of the same rank, ax and xx, is defined as 



(ax ■ Xx) = (-irV2A + l[a A x X a]oo- (A.10) 

It is also convenient to formally introduce the contravariant components a^ of a 
tensor ax by relation 

< = «V (A.ll) 
The scalar product can then be written as 

(«a • Xx) = 5^(-l) M ttA M XA- M = Yl = Yl a ^x' 1 - (A.12) 

The square of a tensor ax expressed in terms of the variables g-a/j and bx^ reads 

A A 

(ax-ax)=J2 a l, + J2 b ^- ( A - 13 ) 

fi=0 fi=l 

This shows that oa's and 6a 's can play the role of the Cartesian coordinates in the 
(2A + l)-dimensional Euclidean space, whose volume element is 

A A 

dQ = Hdax^Udbx,. (A.14) 

fj,=0 fi=l 

Let this Euclidean space be infinite, i.e. — oo < <2a m , bx^ < oo for all /i's. The scalar 
product in the Hilbert space of functions of the Cartesian coordinates a AM , bx^ is 
defined in the usual way by 

($|$) = J **(ax,bxMax,bx)dn. (A.15) 

The integral is over the whole space. 
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Apart from the spherical tensors, the bitensors or matrices with tensor indices also 
appear in the collective models. The spherical bitensor is defined in a similar way as 
the tensor, namely 

M'^ >KU , = Y^ D lV(^^)DTA^^)Mx^ (A.16) 

where M\^ KV and M x ^, , are the components of the bitensor Mx, K in the coordinate 
systems U and U', respectively. The bitensor M\ tK can always be decomposed into a 
number of tensors T^\f of ranks L, where |A — k\ < L < A + k, as follows: 

M Xfi}KU = J2(^\ LM ) T lm- (A-17) 

LM 



If a matrix of the bitensor M\ \ is symmetric the summation in (|A.17I) is restricted to 
even L's. The scalar product of a bitensor M\, K and a tensor ot K is defined similarly to 

(M ViK • a K ) = ^(-l)"M AftK _,a w = ^ M x ^a KV (A.18) 

V V 



The differential operator: 

^ = ^ = ~ ih ^k (A - 20) 

are Hermitian in and play the role of the momenta canonically conjugate to the 
coordinates since they fulfill the following commutation relations 

[ax^Pxu] = ihStxv [bx^Pxu] = [ax^qxA = [bx^QXu] = iMp,. (A. 21) 
The covariant tensor operator of momentum is defined as 

^ = ^(Pa, + ^) = -ih-^- = -iftJ^L. (A.22) 



The relations inverse to (IA.7I) have been used here to express the derivatives with respect 
to otXfi through the derivatives with respect to oa m and &a^- The momentum tensor 
Hermitian adjoint to ir^ is 

*l = ^^-^) = -y^-. (A.23) 

It plays the role of the momentum tensor canonically conjugate to the coordinate tensor 
ax since it fulfills the following commutation relation: 

= ihS ^- ( A - 24 ) 

% Differential operators acting in the Hilbert space of functions of variables cl\^ , are denoted with 
the hat throughout the paper. 
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It follows from (IA.1I ) that the Hermitian adjoint momentum operator ( 1A.23I ) has the 
transformation properties of the contravariant and not the covariant tensor in spite of 
the convention of notation of ( 1 A. 1 lh . According to the convention (1A.11I) 

d_ 

Obviously, the momentum tensor is of the same parity as the conjugate coordinate 
tensor. It is easy to check that the differential operator 



*a" = >»— = -*U (A.25) 



obeys the following commutation relations: 

[L[ X J, a XfM ] = v / A(A + l)(A / ulz/|A / u + v)a^ +u (A.27) 

and 

[L[ x J,7r Xtl ] = v/A(A + l)(A/ilz/|A/i + v)n XfM+u . (A.28) 

This means that L^j is the angular momentum operator acting in the space of functions 
of a Xfi (cf [34]). The angular momentum (1A.26I) is always an axial vector irrespective 
of the parity of a x and n x . Equations (IA.26|) . (IA.27I) and (1A.28|) lead to the obvious 
commutation relations for the components of the angular momentum 

MiUS?] = -hy/2Q.vLlE,\lK)LW. (A.29) 



Appendix A. 2. Isotropic tensor fields 

In continuum mechanics there appears the notion of an isotropic function of variables 
which are the components of a tensor (cf [120]). Roughly speaking, the function 
is isotropic if no parameters with tensor properties (e.g. material tensors) enter 
its definition. The isotropic functions which are tensors themselves, are called the 
isotropic tensor fields. The isotropic tensor field T LM (a x ) of rank L obeys the following 
transformation rule: 

L 

T' LMI {a x ,)=T LM ,{a' Xll )= ^ D L ^ M ,^^^)T LM (a x ,) (A.30) 

M=-L 

where Tlm and a Xfl are the components of tensors Tl and a x in the original system U 
whereas T' LM , and a' X/1 are the respective components in the rotated system U'. Equation 
OA. 301) states that to transform the field Tlm(&\) it is sufficient to transform its tensor 
argument a x ^. It would not be so for a tensor field which is not isotropic. 

The question arises to what extent the form of functional dependence of the isotropic 
tensor field T LM (a x ) on the tensor a x is determined solely by their tensor properties. To 
answer this question we present below a scheme of constructing an arbitrary spherical 
tensor as a function of another spherical tensor. It turns out that for a given rank 
A there exists a finite complete system of i x irreducible (i.e. none of them can be 
expressed rationally and integrally in terms of the others) elementary tensors eu m (a x ), 
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i = of different ranks /. The highest components (m = /) of the elementary 

tensors are called the elementary factors [40]. The elementary tensors are constructed 
by coupling successively n tensors a x to different Ts, I < n\: 

eu m (ax) = e$(a x [c]) = \ a x x ... x a x ] fl, (A.31) 

n 

where the symbol [c], redundant in most cases, specifies the coupling scheme. Tensor 
a\ = e^ 1 ' itself is an elementary tensor. The scalar product of two tensors a x , 

(2) 

(a x ■ a x ) = e as well as all other independent scalars are the elementary factors. 
Completeness of the system of elementary factors does not mean that relations between 
them must not exist. On the contrary, some relations between them, called syzygies, do 
exist. They take the form: 

S j (e lu ,...,e ix w) = Q (A.32) 

with j = l,...,j x , where Sj are some rational integral functions. The number of 
independent syzygies is finite and depends on the rank A. It turns out that to obtain a 
tensor of a given rank L, it is sufficient to align (i.e. to couple to the maximal angular 
momentum) the elementary tensors. The alignment of the elementary tensors means 
the multiplication of the elementary factors (the relevant Clebsch-Gordan coefficients 
are all equal to 1). Because of the existence of the syzygies ( 1A.32I) . some of the aligned 
tensors 

[e{ n) (a A [c]) x e^'\a x [c']) x . . ]l=i+i>+... 

can be expressed rationally and integrally in terms of the others. Using the syzygies 
the finite number k X L of independent aligned tensors r kLM (e) with k = 1, . . . , k XL can 
be found. These independent aligned tensors are called fundamental tensors. This is 
because an arbitrary tensor field of rank L, being the isotropic function of tensor a x , 
can always be written in the form: 

T L m{u\) = ^ fk(so)T kLM (e) (A. 33) 

k=l 

where argument Eq of /& stands for all independent scalars. Scalar functions can be, in 
general, arbitrary. However, the form of some of fk can be restricted. For instance, some 
scalar Eq cannot appear in it in too high power. This is because a syzygy can make the 
expression (£ioo) n (£i'zz) n/ f° r some i, i' , n and nf dependent on other elementary factors. 
When a x has the positive parity, an arbitrary field of the form ( 1A.33I ) has obviously 
the positive parity as well. However, when a x has the negative parity, the field Tl has, 
in general, no definite parity. To decompose it into the parts of the positive and the 
negative parity, the parities of the scalars and the fundamental tensors should be taken 
into account. This can easily be done by noticing that the parity of an elementary tensor 
e<T\[c]) is (-1)». 

A justification of the procedure presented above can be traced back to the theory 
of covariants of algebraic forms given by Dickson |121| . which, however, uses different 
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notions than these used in the present paper. Let us sketch the main points of this 
theory (cf |122l l37j). Let the binary (in variables x\, x 2 ) 2A-ic (of order 2A) algebraic 
form be 

F<z\{xx, x 2 ; a x ~x, • • • , aoa) = E \ ax^x^x^ M (A.34) 

„=_A V A " » / 



with a set of coefficients (acx^, A* = — \ . . . , A). Changing in (IA.34I) the variables by a 
linear non-singular transformation 



2 

Xi 

k=l 

where i — 1, 2, one obtains 



(A.35) 



^(^1,^2;^) = G^a^/i, 2/25 «a) = V] . Xa m («a)?/i +/ '?/2 " (A.36) 

m= -aV A -^/ 

with a new set of coefficients (xxp, I 1 — ~ \ • • • , A). The 2Z-ic form 



1 / \ x / 2 

2 > 



^r ) (x 1) x 2 ;a A )= E ( * m ) h%}(ax)x[ +m x l ->" 



m=—l 



(A.37) 



such that 



H^{yi,y 2 ;xx) = (det{A ik )) w (x u x 2 ; a A ) (A.38) 

is called a (homogeneous) covariant with the weight u> of F 2 \. Coefficients are 
homogeneous polynomials of order n (called the degree of the covariant) such that 
Xn = l + w. For I = the form Hq = is called an invariant of F 2 x- The polynomial 
in front of the highest power of x\ in (IA.37I) . h[] \ is called a semi-invariant of F 2 \. 
It can be shown in the theory of covariants of algebraic forms that any covariant can 
be expressed rationally and integrally in terms of a finite, irreducibly complete set of 
covariants, which can be related rationally and integrally to each other by a finite system 
of independent syzygies (Gordan-Hilbert Finiteness Theorem, see [12 2]). The number 
of the basis covariants for a few lowest A's is listed in [37] . 

It remains to answer the question of what the covariants of algebraic forms have in 
common with the tensor fields being functions of tensors. It turns out that the semi- 
invariant hii(o£\) forms the component of the highest projection of a tensor field of rank 
/ depending on the tensor of rank A (the highest weight of an irreducible representation of 
SO(3) embedded in an irreducible representation of SU(2A + 1)). The other polynomials 
him ( a >) m AA.37I ) are the remaining components of the same tensor. A proof of this 
statement is demonstrated briefly in reference |123| . Constructive proofs of the form 
(1A.33I) can be carried through by explicitly constructing a basis in the space of functions 
of a\ and showing that it has the structure as in (CP3l . In the case of A = 1 it is well 
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known from textbooks of quantum mechanics that an arbitrary tensor field of rank L 
which is a function of the position vector r takes the form: 

T LM {r) = f L (r)Y LM (~) (A.39) 



where r 2 = r r and Ylm is the corresponding spherical harmonic. From (IA.39I) it is seen 



immediately that in the case of A = 1 there are the two elementary tensors [i\ = 2), 
namely the vector r and the scalar r 2 , and no syzygy (ji = 0). Moreover, for every 
given L there is only one fundamental tensor (k\L = 1) of the form [34J: 



r LM (r) = [ rxrx . . . x r] LM cx r L Y LM (-). (A.40) 
v r 

L 

Constructions of the oscillator bases in the cases of A = 2 (see [42] and references quoted 
therein) and A = 3 |123j have been also performed confirming the form ( 1A.33I) . 

Appendix B. Properties of the ATDHFB inertial functions in the intrinsic 
frame 

Appendix B.l. (In) dependence on the Euler angles. 

In this Appendix the three Euler angles (if, ip) will be compactly denoted as cp or ipj, 
j=l, 2, 3. The BCS-type states depending on the deformation and the Euler angles if 
are constructed as £/(<£>) |$(q)) where |$(<?)) is the BCS state in the intrinsic frame and 
q stands for any of the pairs of the deformation variables: (a ,a 2 ), (^0,^2) or (/3,7). 
U((p) is the rotation operator in the Fock space. The argument (f will sometimes be 
dropped for the sake of brevity. The density matrix TZo(q, <p) can be written as 

H Q (q^)=Wl Q (q)U + (B.l) 




u = : . (b.2) 



where 



u is a matrix of U((fi) in a fixed basis of the one-particle space. It fulfills the relation 

V 

Equation (IB.lft can be derived as follows (the (11) block of the TZq is taken as an 
example): 

ffi\q,ip))mn = (U${q)\d+d m \U<S>{q)} = {<$>(q)\U + d+UU + d m U\<$>(q)) = 

= Y,<, u ^{q)\d^d k \<$>{q)) = {uKf{q)u+) mn . (B.4) 

jk 

Now we show that the ATDHFB inertial functions do not depend on the Euler 
angles. To calculate the vibrational functions B qiQk one needs the derivatives of 1Zo(q, V 5 ) 
with respect to the deformation variables. It is easy to show that the derivatives 
can be written analogously to (jB.lj) . i.e. d qj 7Zo(q, (p) = lAd qj 1Zo(q)U + . The same 
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is true for TZ^(q,cp) in view of (14.81) and the properties of the induced self-consistent 
Hamiltonian W. Consequently, all vibrational functions (which are proportional to the 
trace Tr 2 d{Ri [d qj lZo, T^-o])) do not depend on the Euler angles. 

The rotational functions deserve a detailed treatment. Derivatives of 7Z(q, <p) with 
respect to the Euler angles (pj are equal 

VCfavO = {d Vj U<S>{q)\d+d m \U®{q)) + {U<$>{q)\d+d m \d^U<$>{q)) = (B.5) 
= (U-'d^U^lU+d+drnUMq)) + (Q(q)\U + ($d m U\U- 1 dv.U$(q)) = (B.6) 
= {un^{q,ip)u + ) mn (B.7) 

where 

(^Si)«m(9,v) = (HqMd^U-^U^q)) +h.c. (B.8) 
The operator U~ x d v XJ can be expressed in terms of the generators of the rotation group 

tr 1 ^ = ^p fci ( v >)y fc (b.9) 

it 

where P((p) is the universal (i.e. the same for all representations of the rotation group) 
matrix 

' — sin?? cos ip, smtp, ' 
P — sin-^sin^, cos?/', (B.10) 
^ cos-#, 0, 1 j 

and Yfc = —iJ^/H, where are components of the angular momentum operator in the 
Fock space. Finally 

^oifovO = E P ^(v)^ n (9)« + (B.11) 

where the matrices 

(J- fc n (q)) mn = -ift-^^ld+^Jfcl^g)) +h.c. (B.12) 

do not depend on the Euler angles. Of course, all the blocks of the lZoj(q,(p) matrix 
exhibit the same dependence on the Euler angles. Moreover, the matrices TZ\ behave in 
an analogous way, that is 

TZ[ = (pi^2P r i{<p)UZr{q)U + (B.13) 

r 

where Z r (q) do not depend on the Euler angles. 

Taking into account that the angular velocity in the body fixed frame is given by 

^j = ^2Pjk<Pk (B.14) 

k 

the rotational classical Hamiltonian can be written as 

#cl,rot = ^ / ] Ikj^hUj (B.15) 
kj 

with 

/«=*W*RJ) (B.18) 
which obviously do not depend on the Euler angles. 
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Appendix B.2. Block structure of the inertia tensor in the intrinsic frame 

Taking into account the D 2 symmetry of the | $(<?)) state, one can prove the ATDHFB 
inertial functions give the classical kinetic energy which has the form identical with that 

of n 

-Hcl.kin = #cl,kin,vib + #cl,rot (B-17) 

or, more precisely, that the vibration-rotational terms are absent and, moreover, that 
hj — $kjlk ■ (B.18) 

The sketch of the proof is as follows. The operators = exp(— mJ k /K) are the 
generators of the group D 2 and invariance of the \&(q)) with respect to the D 2 group 
means that Sk\<&(q)) = |$(q)), which in the density matrix language can be rewritten 
as S^R-qS^ = TZq. The connection between Sf. and S k is analogous to that between 
U and U (cf ( 1B.ll) to ( 1B.3I) ). For the angular momentum in the Fock space one has 
S^JjSk = Jj for k = j and S^JjSk = — Jj otherwise. In consequence, the matrices Tj 
and Z\ have the same property, i.e. S^TjS^ = The formula (1B.18I) follows from 

the definition (1B.16I ) of Ijy and the properties of the trace. In a similar way one can 
verify that the mixed vibration-rotational parameters B qiipj vanish. 

If H c \ rot ( 1B.15I) is treated as an expression for the metric tensor in the space of 



rotation matrices and (|B.18l) is valid, the corresponding Laplace-Beltrami operator has 
the form 

( B -!9) 

k lk 

where L k are the generators of the right regular representation of the SO (3) group. They 
are proportional to the collective angular momentum operators in the intrinsic frame, 
cf (121511 . 
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